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OIL  SPILLS  IN  THE  ARCTIC 

1.  INTRODUCTION 

With  the  discovery  of  large  oil  reserves  in  the  Arctic, 
there  naturally  arises  questions  about  the  behavior  of  oil  spil- 
led, as  it  may  change  the  fragile  Arctic  environment.  To  date, 
there  have  been  few  oil  spills  in  the  Arctic.  Nevertheless,  one 
expects  that  tankers,  pipelines  and  drilling  will  all  contribute 
to  the  oil  spillage,  much  as  happens  in  warmer  waters.  However, 
as  there  is  little  in  the  way  of  a historical  record  of  oil  spil- 
led in  the  Arctic,  it  is  not  presently  possible  to  determine  which 
of  these  sources  is  the  most  significant. 

The  environmental  concerns  aroused  by  potential  oil  spillage 
are  founded  broadly  on  two  ideas.  First,  the  characteristic 
time  for  biological  degradation  in  the  high  Arctic  is  believed  to 
be  quite  long,  on  the  order  of  ten  years.  This  means  that  if 
spillage  in  an  area  is  more  frequent  than  once  a decade,  we  might 
expect  permanent  changes  in  the  environment.  Second,  crude  oil 
is  generally  a dark  color.  If  spilled  on  ice  or  tundra,  it  will 
change  the  solar  albedo  of  the  area  covered.  Such  a change  in 
solar  heating  is  known  to  drastically  change  tundra.  Of  course, 
there  is  also  the  toxic  effect  of  oil  on  flora  and  fauna  - but 
a discussion  of  such  biological  effects  is  beyond  the  scope  of 
this  report. 

There  are  a number  of  new  features  which  arise  in  Arctic  oil 
spills.  First  is  the  role  of  ice.  Sea  ice  has  a variable  den- 
sity, ranging  from  perhaps  0.70  (gm/cc)  to  0.91  (gin/cc) , the 
latter  value  being  associated  with  old  ice  containing  very  little 
brine.  The  actual  ice  density  depends  on  how  the  ice  formed,  and 
its  subsequent  history. 
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In  contrast,  the  density  of  North  Slope  crude  oil  is  0.89. 
After  aging  two  weeks  in  the  Arctic  Sumner,  this  density  rises 
to  about  0.95^.  Hence  the  oil  can  be  either  heavier  or  lighter 
than  the  sea  ice. 

Suppose  now  that  a large  volume  of  oil  is  spilled  in  a lead. 
Subsequent  .to  the  spill,  the  lead  closes.  On  a hydrostatic  argu- 
ment, one  can  say  that  if  the  oil  is  heavier  than  the  ice,  it  will 
be  forced  under  the  ice  pack.  The  Important  idea  here  is  that  we 
must  consider  both  situations  - oil  under  and  over  the  ice. 

A second  new  feature  is  that  of  surface  roughness.  It  is 
(2) 

well  known  that  sea  ice  is  rough,  and  that  the  roughness 

elements  have  a more  or  less  random  character,  with  "hills"  and 

"valleys".  These  statistical  properties  of  the  Arctic  surface 

are  extremely  variable  - the  rms  roughness  height  can  be  as 

small  as  several  inches  or  as  large  as  several  feet.  Similar 

(2) 

statements  can  be  made  for  tundra  . 

Roughness  is  Important,  because,  if  oil  is  spilled  on  a 
rough  surface,  it  will  run  down  into  the  "valleys"  and  remain  in 
these  "pockets".  See  Figs.  1 and  2.  After  oil  spilled  has  stop- 
ped spreading,  its  area  is  determined  by  the  extent  of  these 
pockets  or  puddles.  Of  course,  the  rough  surface  also  modifies 
the  way  in  which  oil  spreads  over  the  surface. 

These  effects  may  also  be  present  if  oil  spreads  under  the 
ice.  In  that  case,  puddles  form  at  the  high  points  of  the  ice- 
water  Interface. 

A third  feature  is  the  fact  that  crude  oils  are  sticky,  and 
they  stick  to  ice.  This  feature  means  that  ice  traps  oil,  pre- 
venting its  free  release  and  spread  in  the  environment.  Of  course, 
in  the  Arctic  summer,  the  ice  adjacent  to  oil  can  melt,  and  the 
oil  released.  Hence  one  would  expect  that  oil  spilled  in  the 
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FIG.  I GENERAL  ASPECT  OF  ONE  REALIZATION  OF 
THE  PROFILE 


FIG.  2 AVERAGE  LONGITUDINAL  DIMENSION  OF 
A POCKET 
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Arctic  would  be  cyclically  trapped  and  released. 

A fourth  new  feature  has  to  do  with  how  oil  spreads  on  water. 

(3) 

This  is  a topic  which  Is  fairly  well  understood  . The  new 
feature  in  the  Arctic  situation  is  that  the  surface  tension 
spreading  regime ^ is  apparently  absent.  Thus,  Arctic  oil 
spills  do  not  spread  to  a thin  film  (driven  by  surface  tension) 
as  do  their  warm  water  counterparts. 

The  purposes  of  this  report  are  twofold.  First,  we  shall 
give  a description  of  the  short  term  history  of  an  oil  spill  - 
including  the  rate  of  spread,  the  area  covered  by  the  oil,  and 
the  effects  of  leads.  In  each  of  these  areas,  our  knowledge  is 
based  upon  a physically  sound  theory  which  successfully  cor- 
relates carefully  controlled  laboratory  data,  and  what  field 
data  is  available. 

Second,  we  discuss  the  evolution  of  oil,  from  a non-bio- 
logical  point  of  view.  A model  of  oil  aging  is  developed,  and 
its  predictions  compared  with  field  observation.  The  stickiness 
of  oil  on  sea  ice  has  been  measured,  and  the  results  correlated 
with  a theory  which  predicts  how  thick  a layer  of  oil  can  stick 
to  the  ice  in  Arctic  conditions. 

From  these  studies  we  show  that  the  evolution  of  an  oil 
spill  in  the  Arctic  is  a much  slower  phenomena  than  in  temperate 
waters.  The  area  covered,  even  from  a supertanker  spill,  is 
quite  small,  on  the  order  of  some  few  square  kilometers.  Thus 
the  technological  problem  of  combatting  Arctic  oil  spills  is 
much  different  than  that  in  temperate  water:  it  is  not  likely 

that  barriers  will  be  useful  or  needed.  A system  with  rapid 
response  time  is  not  required.  On  the  other  hand,  there  seems 
to  be  no  presently  known  way  to  locate  easily,  much  less  remove, 
oil  trapped  under  oil  in  ice. 
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2.  DESCRIPTION  OF  SPILL 


In  thi-  section  we  give  an  analysis  of  the  events  Immediately 
following  the  release  of  oil  In  the  Arctic.  To  do  this  we  must 
develop  a description  of  the  rough  Ice  surface.  We  shall  con- 
sider the  ice  to  have  an  average  height  <h>,  above  sea  level. 

The  actual  height  of  the  ice  is  then  equal  to 


! 


h(x)  ■ <h>  + h’ (x) 


(1) 


h'(x)  denotes  the  variable  roughness  of  the  Ice  surface. 

Where  h'  Is  negative,  oil  spilled  will  tend  to  accumulate. 

(3) 

The  data  on  h'  are  typically  gotten  from  a laser  pro- 
f Home  ter  mounted  on  an  aircraft  flying  over  ice.  The  signal 
which  results  from  such  a measurement  appears  essentially  ran- 
dom, There  seem  to  be  no  features  which  distinguish  one  area 
of  rough  ice  from  another  (nor  should  there  be,  on  physical 
grounds).  If  no  statistical  distinction  can  be  made  between 
different  areas  of  rough  ice,  then  mathematically  speaking  the 
random  variable,  h’(x)  may  be  considered  a stationary  stochastic 
function. 

To  study  the  pocketing  of  oil  on  the  rough  ice  surface, 
one  is  lead  to  ask  over  what  extent  h'  is  negative,  so  as  to 
determine  how  large  the  pocket  is. 

But,  the  standard  analysis  of  ice  roughness  amounts  to  com- 
puting the  power  spectrum  of  h'.  The  power  spectrum  P is  the 
Fourier  transform  of  the  auto  correlation  function,  ( f is  a dunmy 
variable) 

f 

<h'(x)  h’(x  + t)>  ■ V P(f)cos  2 f df  (2) 

^-00 
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We  can  make  such  a statement  because,  as  no  region  of  the  sur- 
face Is  statistically  distinct  from  other  regions,  the  average, 
<h'(x)  h'(x  + t)>  only  depends  on  T.  Then  P only  depends  on  f. 

Up  to  now,  we  have  been  discussing  how  one  region  of  the 
surface  is  correlated  with  another,  located  at  an  arbitrary 
distance  T away.  However,  the  actual  data  are  now  more  re- 
strictive, as  the  aircraft  flies  straight  while  the  prof ilometer 
records  surface  roughness.  For  our  application  T is  a distance 
along  the  surface  parallel  to  the  flight  path. 


<h'(x)  h' (x  + t)>  where  h*(x) 


^h'  (x)  , h < 0 
.0  elsewhere 


(3) 


This  is  a major  theoretical  problem.  In  this  report,  we  shall 

(4) 

only  summarize  the  mathematical  theory  we  have  ueveloped  . 

A detailed  physical  description  of  the  results  will  be  given. 


It  turns  out  that  to  relate  (3)  to  (2),  additional  informa- 
tion, beyond  P(f)  is  required.  Basically,  we  use  the  fact  that 
h is  normally  distributed,  i.e.  if  6(h)d’n  gives  the  fraction  of 
the  time  h will  be  between  h and  h + dh,  then  6(h)  has  the  form 


6(h) 


1 

exp 

/2  TO. 

n 


-(h  - <h>)" 

20. 


(4) 


2 

Here,  0.  ■ <(h')  >,  the  mean  square  roughness  height.  Figure  (3) 
n 

fh 

compares  i 6(x)dT  against  h for  a suitable  value  of  0^.  The 
y_oo 

solid  line  is  the  theory,  and  the  points  are  field  data  from 
reference  (3). 
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VERIFICATION  FOR  h NORMALLY  DISTRIBUTED 


The  remaining  parameters  which  characterize  the  surface 

are  measures  of  the  horizontal  scales  of  the  surface.  We  have 

(3) 

found  that  the  field  data  fit  a power  spectrum  calculated  by 
assuming 


<h'(x)  h'(x  + t)>  ■ — 5-  cos  (5) 

1 + kT  1 

provided  suitable  choices  of  k,  T and  c ^ are  used.  Figures  (4, 

5 and  6)  show  the  comparison  between  field  data  and  Eq.  (5). 

Physically,  ^57  Is  the  roughness  height.  T is  the  horizontal 

-1/2 

length  of  a given  pocket,  k is  the  distance  in  which  the  cor- 

-1/2 

relation  between  h'(x)  and  h’(x  + T a x + k ) becomes  negligible. 

Suppose  now  that  a surface,  characterized  by  o^,  k and  T is 
specified.  Let  oil  be  released  on  this  surface.  The  oil  spreads, 
and  we  next  proceed  to  estimate  how  fast  the  oil  spreads. 

There  Is  little  doubt  that  the  oil  spreads  due  to  gravity. 

Let  capital  H be  the  thickness  of  the  oil  layer,  and  A the  dif- 
ference In  specific  gravity  between  oil  and  water.  Then  the 
driving  force  Is 


2x(pAgH)Hr  where  A 


C l over  ice 


^A  under  ice 


(6) 


as  r is  the  radius  of  the  oil  pool  at  a time  t. 

After  careful  study  of  all  available  data,  we  have  found 
that  the  retarding  force  is  a frictional  one,  as  the  oil  reaches 
new  pockets.  This  frictional  drag  Is  concentrated  at  the  leading 
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<rh  = 1.016  ft2 

k = 0.002786  ft2 
T = 200  ft 


edge  of  the  oil  slick.  For  radii  less  than  r,  the  oil  flows 
over  filled  pockets  of  oil  with  essentially  no  friction. 

Because  this  pocket  filling  process  is  controlling  the 
rate  at  which  the  oil  spreads,  it  turns  out  that  the  oil  spreading 
rate  is  controlled  by  rate  of  release  of  the  oil.  For  simplicity, 
we  shall  assume  an  average  rate  of  volume  released,  0.  Continuity 
then  requires  that 

*r2H  - Qt  (7) 

The  frictional  retarding  force  is  proportional  to  the  "area" 
of  the  roughness-  reen  by  the  approaching  oil,  (<h^>  - — 


2irr<h^> 


The  pressure  drop  experienced  as  the  oil  flows  into  new  pockets 
is  proportional  to,  and  is  caused  by  the  drag  exerted  by  the  rough 
surface  on  the  oil. 


P 


( — ) 
v dt  1 


2 


where  is  the  velocity  of  the  oil  at  the  leading  edge  of  the 
slick.  Balancing  the  driving  and  retarding  force  gives 


2ff(pAgH)Hr  ~ 2TTr<h’>  p 


v dt  ' 


(8) 
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By  using  Eq.  (7),  we  can  determine  that,  as 


dr  r 
dt  “ t* 


(the  coefficient, 


0.25  is  determined 
r 


empirically) 


<9) 


Fig.  (7)  shows  that  this  simple  equation  correlates  both  the 
laboratory  data^  and  field  data^’^.  In  working  out  calcu- 
lations of  how  fast  oil  spreads  over  or  under  ice,  Eq.  (9)  should 
be  used. 

The  laboratory  data  were  gotten  by  building  several  random 
rough  surfaces  whose  characteristics  could  be  fitted  empirically 
by  Eqs.  (4  and  5).  Experiments  were  done  both  for  oil  spreading 
under  and  over  the  ice  surface^. 

Oil  spreading  in  a lead  filled  with  water  is  easily  treated 

as  are  other  spills  over  water That  is,  there  are  several 

phases  of  oil  spreading  - inertia,  viscous  and  finally  surface 

(8) 

tension.  There  are  standard  formulae  which  give  the  rates 
of  spread.  However,  there  is  one  important  modification  which 
arises  in  the  Arctic;  the  surface  tension  phase  of  oil  spread 
is  absent,  because  the  net  spreading  coefficient  is  essentially 
zero.  This  means  that  the  motion  ceases  when  the  oil  reaches 
a certain  thickness,  on  the  order  of  lcm^.  Practically  speak- 
ing,  then,  Eq.  (9)  and  those  referenced'  ' provide  a method  of 
working  out  the  spreading  rate  of  oil  spilled  in  the  Arctic. 

The  spreading  stops  when  either  (a)  the  oil  thickness  on 
water  reaches  about  1cm,  or  (b)  if  the  oil  fills  all  the  pockets 
of  the  ice  field.  In  the  first  case,  the  area  of  water  covered 
by  oil  is  simply  proportional  to  the  volume  released: 


Area  ■ Volume/thickness 


(10) 
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FIG.  7 DIMENSIONLESS  RADIUS  vs  TIME  FOR  LABORATORY  AND  FIELD  DATA 
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The  situation  is  more  complex  in  the  second  case.  If  the 

(A) 

rough  Ice  data  Is  fitted  with  Eq.  (5),  then  Patureau  shows 
that  the  surface  area  of  an  average  pocket,  Smp,  is 


Smp 


(ID 


and  the  average  volume  of  a pocket  Is 


Vmp  ■ — j (2ah)1/2  (12) 

IT 

Now  the  area  that  an  oil  slick  will  occupy  depends  simply  on  the 

average  number  of  pockets  which  can  be  filled  with  oil.  Clearly 

the  total  area  covered  must  then  be  proportional  to  the  volume 

(A) 

spilled,  and  inversely  proportional  to  Vmp.  The  exact  result 
is 


<S  > 
max 


Snp  ~jn  ( > 


(13) 


<S  > is  the  average  area  of  a spill  over  or  under  ice. 
max 

For  a random  rough  surface,  Figs,  (8  and  9)  compare  this  result 
with  laboratory  data.  It  is  seen  that  there  is  fairly  good 
agreement  with  this  result. 

Using  this  result,  we  can  estimate  the  area  covered  by  an 

oil  spill  on  ice.  Table  I shows  the  results  for  a spill  of 
5 3 

V ■ 5 x 10  ft.  (about  15,000  tons)  for  three  sets  of  T,  a, 

(A)  n 

gotten  from  Arctic  data  . 
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FIG.  8 DIMENSIONLESS  VOLUME  RELEASED  vs  SURFACE 

AREA  FOR  SURFACE  1 
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TABLE  I 


T 

0h 

Smp 

<S  > 
max 

115  ft. 

2,295  ft.2 

1.32  x 104  ft.2 

7.29  x 105 

ft. 

200  ft. 

1.016  ft.2 

4 x 104  ft.2 

1.24  x 106 

ft. 

540  ft. 

0.203  ft.2 

4 x 104  ft.2 

2.78  x 106 

ft. 

If  a super  tanker  full  of  oil,  150,000  tons,  were  spilled, 
then  for  ice  described  in  Table  I,  one  gets  areas  ranging  from 
1/4  to  1 square  mile  as  the  area  covered.  Such  small  areas 
should  be  contrasted  with  the  Torrey  Canyon  experience  - in  that 
case,  the  oil  spread  to  some  300  square  miles  in  area.  On 
Arctic  ice,  in  contrast,  the  present  study  predicts  a much 
smaller  area. 

The  calculations  of  <S  > in  Table  I are  based  on  one 

max 

dimensional  spectra,  and  the  assumption  that  the  statistics  are 
the  same  in  all  directions  on  the  ice  surface.  This  is  a crude 
assumption,  which  was  made  because  no  two  dimensional  spectra 
of  ice  field  roughness  are  currently  available.  However,  using 
a variety  of  assumptions  of  the  asymmetry  of  ice  roughness 
elements  does  not  change  the  result  by  much  more  than  a factor 
of  three. 

Although  there  is  very  little  field  data  on  the  character- 

(3) 

lstics  of  the  under  ice  roughness  , it  is  approximately  true 
that  the  ice  is  in  hydrostatic  equilibrium.  If  we  assume  that 
this  is  exactly  true,  then  for  the  under  ice  surface,  T is  the 
same  as  that  of  the  over  ice  roughness,  and  the  roughness  height, 
is  about  seven  times  greater;  as  1/7  of  the  ice  is  above 
the  mean  water  level.  Thus,  the  mean  pocket  volume,  Vmp,  is 
seven  times  greater  for  the  under  ice  surface.  (See  ref.  2) 
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If  we  consider  two  cases,  one  in  which  a volume  V of  oil 

is  spilled  over  the  ice,  and  another  where  the  same  volume  is 

spilled  under  the  ice,  then  the  above  argument  suggests  that  the 

area  <S  > of  the  oil  under  the  ice  is  1/7  that  of  the  area 
max 

over  the  ice. 

The  hydrostatic  assumption  is  too  crude  to  permit  the  ac- 
curate estimate  of  this  area  ratio.  But  it  is  quite  clear  that 
the  Vmp,  mean  pocket  volume,  under  the  ice  must  be  much  greater 
than  over  the  ice. 

All  the  results  of  this  section  may  be  summarized  as  follows: 
At  the  end  of  the  spreading  phase  of  an  Arctic  spill,  the  oil 
will  be  trapped  either  in  a lead  or  on  the  ice  surface;  and  the 
area  will  be  much  less  than  what  one  would  expect  in  temperate 
climates.  Practically  speaking,  we  expect  that  a super  tanker 
spill  will  occupy  at  most  about  a square  mile  in  area. 

3.  EVOLUTION  OF  OIL  IN  THE  ARCTIC 

In  this  section  we  consider  what  happens  after  the  spilled 
oil  has  ceased  to  spread.  Then  the  oil  will  either  be  in  pockets 
above  or  below  the  ice  surface,  or  in  a lead,  floating  on  sea 
water.  From  a mechanical  point  of  view,  we  are  interested  in 
whether  the  oil  may  migrate,  over  a long  period  of  time,  from 
the  origin  of  the  spill.  There  are  really  only  two  potential 
causes  for  such  migration.  First,  there  is  the  slow  increase 
in  oil  density  with  time  - we  refer  to  this  process  as  aging. 

In  about  two  weeks,  in  the  Arctic  summer,  ^ the  density  of 
Forth  Slope  crude  increases  to  0.95,  which  is  heavier  than 
virtually  all  Arctic  ice. 

In  the  Arctic  summer,  the  dark  oil  on  the  surface  of  the 
ice  absorbs  enough  heat  to  float  in  a pool  of  melt  water.  If 
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a lead  now  opens,  Che  oil  will  flow  into  it.  If  the  lead  then 
closes,  the  oil  will  flow  either  under  or  over  the  ice  depend- 
ing on  whether  it  is  lighter  or  heavier  than  the  surrounding 
ice.  In  this  sequence  of  events,  the  oil  density  determines 
whether  the  oil  will  flow  under  or  over  the  ice. 

The  second  major  factor  in  migration  of  oil  is  the  question 
of  whether  oil  sticks  to  ice  or  not.  We  have  done  a separate 
experiment  to  determine  the  stickiness  of  North  Slope  crude  and 
diesel  oil  to  sea  ice,  and  find  that  substantial  amounts  of  oil 
stick  to  ice.  In  particular,  this  stickiness,  discussed  below, 
apparently  precludes  the  dispersion  of  oil  under  the  ice  by  the 
action  of  currents. 

Because  of  these  somewhat  unexpected  results,  the  evolution 
of  oil  spilled  in  the  Arctic  will  generally  consist  of  an  evapora 
tion  of  the  lighter  fractions,  causing  an  Increase  in  density. 
Second,  if  the  oil  is  located  under  the  ice,  it  will  stick  to 
the  ice.  In  the  winter  season, oil  located  over  the  ice  will 
stick  to  the  ice  also.  In  summer  , the  oil  located  over 
the  ice  does  not  stick,  but  floats  in  a pool  of  melt  water. 

From  this  discussion,  it  is  apparent  that  the  evolution 
of  spill  oil  hinges  on  an  analysis  of  how  the  density  of  the 
oil  changes  with  time,  and  how  the  oil  sticks  to  the  ice  surface. 

A.  The  Evolution  of  the  Boiling  Point  Distribution 

The  purpose  of  this  section  is  to  develop  equations 
which  describe  the  boiling  point  distribution  of  crude  oil 
as  it  ages  on  ice  or  water.  We  assume,  firstly,  that  at  all 
times  the  oil  is  well  mixed.  That  is,  the  concentration  of 
the  ith  component  of  oil  is  always  uniform  throughout  the 
oil.  The  argument  for  this  assumption  goes  as  follows: 


UI  rjm.flywEv  mn^ii  >i!_f ,1'^r^TW 
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If  the  oil  spreads  to  a thin  (1  era)  layer  on  water, 
lntermolecular  diffusion  will  smear  out  any  gradients  of 
concentration  in  a short  tine  (1  day).  On  the  other  hand, 
if  the  oil  is  pocketed  in  deep  pools,  the  action  of  wind 
will  be  to  make  waves  on  the  oil  surface  which  will  mix  the 
oil.  Even  in  the  Arctic  winter,  when  the  oil  on  the  ice 
surface  is  frozen  solid,  we  believe  this  will  be  a good 
assumption  because  the  time  constant  for  evaporation  in- 
creases exponentially  with  decreasing  temperature,  whereas 
the  mixing  time  for  oil  is  at  most  the  time  for  character- 
istic changes  in  the  ice  field,  i.e.  less  than  a month. 

The  assumption  is  that  this  mixing  process  will 
always  have  a time  constant  which  is  less  than  the  time 
constant  for  evaporation. 


Second,  we  use  a simple  turbulent  diffusion  model, 
based  upon  a neutrally  stratified  air  flow  over  a rough 
flat  plate.  There  is  no  question  that  there  will  be  de- 
partures from  neutral  stability  in  the  air  flow  over  the 
oil  or  the  Arctic  Ocean.  However,  these  departures  will 
be  unimportant  near  the  rough  surface,  as  there  the  rough- 
ness elements  will  dominate  the  mixing  process,  provided 

(9) 

there  is  an  appreciable  wind  velocity  . Under  this  as- 
sumption, the  diffuslvlty  at  sea  level  is 


ku*z 


where  z is  the  roughness  height,  u is  the  friction  velocity 
and  k is  Von  Harman' s constant,  having  a value  of  4/10. 


To  compute  a flux  we  need  to  know  the  gradient  of  the 
concentration  of  the  oil  vapor  in  the  air.  We  argue  that, 
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due  to  the  low  rate  of  evaporation,  there  will  be  essentially 
zero  vapor  concentration  for  heights  greater  than  z.  On  the 
other  hand,  due  to  the  strong  mixing  Induced  by  the  rough- 
ness elements,  we  must  expect  finite  concentrations  for 
heights  less  than  the  roughness  height.  Hence,  the  length 
scale  for  the  concentration  gradient  is  z.  The  flux  is  pro- 
portional to  the  concentration  times 


Third,  following  Blockker^^,  we  assume  the  oil,  which 
is  a mixture  of  a large  number  of  organic  components,  is 
considered  an  ideal  thermodynamic  mixture.  Although  this 
assumption  is  certainly  untrue  in  detail,  we  expect  that 
the  discrepancies  caused  by  non-ideal  behavior  will  not 
cause  any  quantitative  changes  in  the  results.  This  assump  - 
tion implies  Rauolt's  law  is  valid. 

In  the  same  vein,  the  oil  vapor  of  the  ith  species  is 
considered  a perfect  gas.  The  vapor  pressure  of  the  ith  com- 
ponent over  its  liquid  obeys  a Clapeyron  equation.  These  are 
excellent  assumptions  because  the  temperatures  of  the  vapors 
are  all  low  compared  to  the  critical  temperature  for  the 
substances  involved,  and  the  gases  are  all  dilute. 

Finally,  let  us  show  that  evaporation  into  the  air 
dominates  dissolution  of  oil  into  the  water.  The  mass  flux 
per  unit  area  of  the  ith  fraction  of  oil  to  the  air  is 

* 
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where  is  the  density  of  the  ith  component  of  oil  vapor 
at  the  oil  surface.  The  flux  of  the  ith  component  to  the 


water  is  about 


where  is  the  solubility  of  the  ith  component  in  sea 
water. 

* * 

Now  u^  is  much  smaller  than  u . To  see  this , consider 
the  case  where  the  oil  lies  at  the  water-air  interface.  At 
the  interface,  the  turbulent  stress  is  continuous.  Hence 


* * /p 

u ■ u / Kalr 


Using  these  results,  the  ratio  of  evaporation  to  air  to  flux 
to  water  is 


o -1  -4 

Typical  values  at  0 C are  c^  * 10  p ^ “ 10  and 

pH20/p«ir  ' 1q3- 

Hence  the  ratio  is  on  the  order  of  30.  He  shall  ignore  dis- 
solution of  oil  into  water. 


3 


Consider  now  a mass  of  oil,  containing  1 fractions, 

each  vith  a molecular  weight  M^.  The  mass  of  each  fraction 
is  m.  The  mole  fraction,  depends  on  the  number  of  moles 
In  the  1th  fraction,  n^,  as 


"1  “1 

T " 


The  mass  of  the  1th  fraction  Is 


"Vi  ’ -i 


Since  the  volume  occupied  by  each  fraction,  is  related 
to  the  density  of  each  component  as 


mi  ■ pivi 


It  can  be  seen  that  the  bulk  density  is 


k m 1 m EMiXl^Pi 

p m ^Mixi 


Raoult's  lawv  ' states  that  the  partial  pressure  of 
the  ith  component  is  proportional  to  Its  saturated  vapor 
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pressure  and  the  mole  fraction  x^ 


P - x P 
ri  xlri 


The  saturated  vapor  pressure  Is  determined  by  Clapeyron's 
equation 


■°  . i)i 

7 • t’J 


This  equation  relates  the  1th  component  to  Its  boiling  tem- 
perature Tg.  At  Tg,  the  1th  component  has  a saturated 
vapor  pressure  of  P0  - 1 atmosphere.  Here  q is  the  latent 
heat  of  vaporization  and  R,  the  universal  gas  constant. 

Next,  we  relate  the  partial  pressure  Pi  to  the  density 
of  the  vapor  about  the  oil,  c^  assuming  the  vapor  to  be  a 
perfect  gas. 


ptMi 


Using  (14)  and  (18)  in  combination  we  determine  the 
rate  of  evaporation  of  the  ith  component  to  be,  for  an 
exposed  area  A, 


>v 


fe  (iH-i)  =j 
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The  solution  of  (19)  gives  n.  (t;T  ) and  N(t).  By  equa- 
tlons  (14)  and  (IS) , this  is  equivalent  to  knowing  M(t)  nd 
Mi(t,TB).  The  evolution  of  the  boiling  point  distribution 
is  normally  represented  In  graphical  form  with  one  axis 

being  the  fraction  having  boiling  temperatures  higher  than 

* 

Tfi,  that  is 


and  the  other  axis  being  T_.  (t*  is  a dummy  variable) 

D O 

B.  The  Evolution  of  Density  in  an  Oil  with  a 
Continuous  Distribution  of  Boiling  Points 

The  purpose  of  this  section  is  to  extend  the  previous 
results  to  an  oil  of  infinitely  many  components,  each  with 
a mole  fraction  dx.  To  do  this,  we  must  associate  with 
each  boiling  point  a molecular  weight,  M(T„).  Of  course, 

D 

there  is  not  a one  to  one  relationship  between  molecular 

structure  of  each  component  of  the  oil.  But  it  is  not 

realistic  to  identify  each  of  the  large  number  of  components 

in  a typical  crude  oil.  Figure  10  shows  the  relationship 

(12) 

M(T_)  for  various  classes  of  hydrocarbons  . The  solid 
line  shown  is  simple  linear  fit  for  the  lower  boiling  points 


M(Tb)  - 0.42(Tb  - 106*) 


(21) 


(Note  that  Tfi  is  absolute  temperature) 
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Knowing  the  molecular  weight,  we  now  proceed  to  con- 


sider an  oil  with  incremental  mole  fractions,  dx(T_,t). 

B 


The  mass  of  oil  is 


f 

H - N I 


MdX(Tg) 


(22) 


where  the  integration  is  carried  out  over  T_.  This  conven- 

B 


tion  will  apply  below. 


The  boiling  point  distribution  may  be  defined,  using 
Eqs.  (22)  and  (20)  as 


i; 


oo 


Mdx 


b<V*>  ’ ft 


I 


(23) 


Mdx 


As  the  number  of  moles  in  each  species  is  dn  ■ Ndx,  Eq.  (19) 
may  be  written 


fe-f 


B 


lu  *p. 


dx  - - 


RT 


r 


B 


f qM  . T " TB 
exp  ( ( — 


)]dx  (24) 


To  evaluate  Eq.  (24)  we  need  to  know  how  q depends  on  T_. 

B 


Trouton’s  rule  (see  Ref.  (10))  is  that  is  approximately 

KiB 

a constant  - its  value  is  near  5.0. 
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Given  b(Tg,  t ■ 0)  and  N(t  ■ 0)  Eq.  (23)  nay  be  used 
to  compute  x(Tfil  t - 0).  Eq.  (24)  plus  the  constant, 


dx  - 1,  x > 0, 


then  specify  how  N(t)  and  (T„,  t)  change  with  time. 

15 

determines  the  evolution  of  b(T_,  t) . 

D 

The  bulk  density  is  related  to  x by 


r1® 

I Mdx 

'T,  Pi*V 


If  all  fractions  which  have  boiling  points  less 

than  T'  have  been  evaporated,  then 


fT 

t:  pi 


Recognizing  that  the  density  changes  are  always  small,  to 
a good  approximation,  in  the  range  50°C  < T_  < 500°C, 


p(T’)  - p(Tg  - 273°)  + < ) (Tj  - Tj  (26) 

B 

where  4iv  ■ 2.8  x 10  ^ , from  Figure  11. 

B ccK 
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For  Kuwait  crude  P(T*  ■ 273°)  ■ 0.869.  In  figure  1 of 

reference  (12)  we  find  that,  for  Kuwait  crude  b(T_)  is 

15 

approximately  a linear  function  of  T ; rising  from  Oat 

D 

50° C to  1 at  500° C. 

b(V  “ (TB  ‘ 323°K>/450°K  (27) 

By  differentiating  Eq.  (23)  and  using  this  result,  we 
determine  that 

(Sf-  " " t ’ A “ const  (28) 

B 

With  this  result  we  may  invert  Eq.  (27)  to  find  the  density, 
Pj(To)t  of  the  fraction  of  oil  which  boils  at  T_.  The  re- 

1 b D 

suit,  obtained  by  using  Eqs.  (24),  (25),  (26),  and  keeping 
only  first  order  changes  in  Ap  is 

Ap..  Ap(T*) 

pi(V  “ pi(Ti>  + (Ap/AT)  <T  - V ; fir  - — xf-  (29) 

Provided  p(T')  and  b(T_)  are  known  for  a given  crude 

D O 

oil,  p,  (T_)  may  be  determined  in  a manner  analogous  to  the 
1 15 

derivation  of  Eq.  (27). 

At  this  point,  we  see  that  if  b(T_,  t - 0),  N(t  - 0) 

15 

and  p(Tl,  t ■ 0)  are  given  we  may  perform  the  following 
steps  to  obtain  a solution  to  the  problem.  First,  p^(Tg) 
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ii  gotten  from  (T^,  t - 0)  and  b(TB>  t - 0) . Second, 

Eqs.  (24)  and  (25)  are  solved  to  yield  N(t)  and  x(T_,t). 

D 

Third,  b(Tg,t)  is  evaluated  using  Eq.  (23).  Fourth,  p(t) 
is  determined  using  Eq.  (29).  In  this  manner,  the  evolu- 
tion of  the  oil  is  specified. 


Calculations 

For  the  purpose  of  calculation,  it  is  simplest  to  de- 
fine the  characteristic  time,  t,  as 


t<tb> 


x"‘pba 

RT 


exp 


T - T 
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RT„  v T 


(30) 


Figure  12  shows  the  behavior  of  the  characteristic  time  for 
various  points  and  various  values  of  T.  In  this  figure  Trouton’ 
rule,  qM/RT.,  * 5.0,  was  used.  This  chart  shows  the  signifi- 
cant  decrease  in  aging  rate  of  volatile  components  as  T,  the 
ambient  temperature,  is  decreased. 

The  first  calculation  we  present  is  for  a two  component 
mixture.  The  lighter  component  has  a boiling  point  of  200°C 
and  the  heavier,  a boiling  point  of  400°C.  For  North  Slope 
crude,  about  25  percent  of  the  oil  has  boiling  points  lower 
than  200°C  and  25  percent  has  boiling  points  greater  than 
400°C.  The  molecular  weights  at  these  boiling  points  are 
M^(200°C)  - 150,  M2(400°C)  = 350  from  Figure  10*  The  den- 
sities are  determined  from  Eq.  (26)  and  p(T'  ■ 273)  ■ .89 

D 


to  be  ■ .94,  p^  ■ 1.0,  at  T 


273. 
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NONDIMENSIONAL  TIME  CONSTANT. rXU  P„A/R  273°K 


w 


Under  these  assumptions  the  equations  of  Section  3 take 
the  simple  form  (ftom  eqs.  14  - 20) 


n^  + n2  ■ N 


1_  \ 
Tl  N 


1_  ^2 

t2  n 


By  dividing  one  equation  by  the  other,  we  obtain 


"l  Z1  , n2  /2 

( nx(t  - 0)  } “ ( n2(t  - 0)  > 


where  the  Initial  values,  n^(t  ■ 0)  and  n2(t  ■ 0)  are  deter- 
mined by  the  assumption  that  both  fractions  have  equal  mass 
at  t ■ 0: 


Vi  " M2n2  • 


which  gives,  for  one  mole  of  mixture  at  time  t ■ 0,' 


nx(t  - 0) 


M~-n*2  ’ n2<fc  ■ 


Mx  + M2 


The  evolution  for  t > 0 is  determined  by  Eq.  (30); 


hence 


» ni<£) 

J nx(t  - 


<n1  + n2<0)) ( 


nl 

n2(0) 


VT2 


By  carrying  out  these  integrals  it  can  be  shown  that  both 
fractions  completely  evaporate  in  the  same  time,  t (see 
Ref.  2) 


max 


t2  n2^ 


0)  + T2  n^(t  - 0) 


and 


t - 


£max  - T1  "l(t>  - T2  n2(t) 


(34) 


Even  though  the  n^  fraction  evaporates  much  more  rapidly 
initially . 

Figure  13  shows  the  boiling  point  distribution 
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The  evolution  of  density, 


1 

P 


(36) 


Is  shown  on  Figures  15  and  14. 

These  results  bear  some  discussion.  On  Figure  13 
notice  that  there  is  a nearly  linear  decay  in  the  light 
fraction,  b^,  on  a time  scale  T^.  This  decay  produces  a 
nearly  linear  increase  in  density  as  shown  in  Figures  15 
and  14.  The  two  component  theory  is  too  crude  to  predict 
the  exact  shape  of  the  density  vs.  time  curve,  as  observed. 

To  compare  these  results  with  observation We  esti- 
mate the  area  the  oil  sample  will  occupy,  and  a typical, 

A 

average  value  of  u . Fairly  good  agreement  is  obtained  be- 
tween theory  and  observation  if  X ■ k/10. 

The  time  scale  for  the  evolution  of  the  density,  t^,  is 
about  10  days  for  T - 273°K  (0°C) , and  about  14  days  for 
T ■ -17°C,  the  temperatures  corresponding  to  summer  and 
winter  conditions.  Much  larger  variations  in  are  pre- 
dicted with  variations  in  wind  speed.  The  prediction  is 
that  during  high  winds,  the  aging  process  accelerates.  This 
is  observed^.  (See  Figure  15).  No  account  is  taken  of  the 
effect  of  snow  cover,  but  when  the  oil  was  covered  with  snow, 
the  aging  process  stopped^.  (See  Figure  15). 

To  account  for  the  variation  of  with  temperature, 
Figure  12  should  be  used.  Using  Figure  13  we  then  can  de- 
termine that,  for  components  boiling  at  200°C,  the  time 
constant  is  increased  from  10  days  to  17  days  when  the 
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FIG.  if 


ambient  temperature  is  decreased  to  -20°C.  From  Figures  IS 
and  14,  it  is  clear  that  is  the  time  to  evaporate  the 
lighter  components  of  the  oil.  Thus,  we  may  say  that  in 
the  Arctic  winter,  the  aging  time  for  oil  is  approximately 
doubled.  Since  it  takes  about  14  days  for  most  of  the  lighter 
fractions  to  evaporate  in  the  summer,  it  take  about  1 month 
for  the  lighter  fractions  to  evaporate  in  the  Arctic  winter. 

We  can  show  that  having  a continuous  distribution  of 
fractions  with  various  boiling  points  does  not  change  the 
picture  outlined  here.  Differentiating  Eq.  (24)  with  respect 
to  Tg  produces  an  equation  similar  to  Eq.  (31).  By  analogy 
with  the  solution  (34)  for  the  two  component  model,  it  can 
be  shown  that  the  solution  to  Eq.  (24)  is 


Jt  T(TB>  % ltdIB  - B [ T<V  % It' 


Since  N dTg  ■ dTg  is  the  moles  in  the  fraction 

B B 

whose  boiling  point  is  T_ , this  result  may  be  written  as 

O 


t - t 


t(Tb)  dn(t , Tg) 


which  can  be  readily  recognized  as  exactly  analogous  to 
Eq.  (34). 

Thus,  we  may  conclude  that  each  fraction  of  oil  decays 
as  described  by  the  simple  two  component  model.  Thus, 
Figure  12  can  be  used  to  predict  the  evaporation  time  of  a 
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component  with  an  arbitrary  boiling  point , T-. 

D 

D.  The  Maximum  Size  of  an  Oil  Spill 

In  this  section,  we  summarize  the  results  of  Wolfe  and 

(13) 

Hoult'  The  detailed  report  of  this  portion  of  the  work 

has  been  written  up  as  a separate  paper,  and  is  included  as 
Appendix  A. 

If  one  assumes  that,  over  a long  period  of  time,  the 
oil,  due  to  density  changes,  migrates  under  the  ice,  then 
the  question  arises,  does  the  oil  stick  to  the  under  ice 
surface?  To  answer  this  question,  an  experiment  was  per- 
formed, in  which  sea  water  was  frozen  by  a one-dimensional 
heat  flux,  thus  simulating  the  Arctic  situation.  Oil  was 
injected  under  the  ice,  and  subsequent  events  were  studied. 
The  heat  flux  was  maintained  after  the  oil  was  injected. 

It  was  found  that,  after  some  time,  the  ice  began  to 
freeze  under  the  oil.  If  the  freezing  process  continued 
long  enough,  an  inclusion  of  oil,  in  ice,  was  formed.  These 
results  directly  imply  that  if  oil  is  spilled  under  the 
Arctic  ice  in  the  fall,  it  will  be  frozen  into  sea  ice 
floes  as  an  inclusion  during  the  winter.  In  this  way,  a 
mixture  of  oil  and  sea  ice  is  generated.  Such  mixtures 
have  been  observed  during  the  Chetabucktoe  Bay  spill. 

Based  upon  this  laboratory  verification,  we  expect  that 
most  oil  spilled  in  the  Arctic  will  end  up  as  in  a mixture 
of  ice  and  snow. 

Another  very  interesting  observation  is  that  North  Slope 
crude,  at  the  temperature  of  freezing  sea  water,  is  a good 
insulator  compared  to  the  sea  ice  itself.  This  was  verified 
by  measuring  the  temperature  profile  in  the  ice  before  and 
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after  the  oil  was  spilled v . This  result  means  that  the 

heat  will  tend  to  flow  around  the  oil  pocket  under  the  ice, 
rather  than  through  it.  Near  the  edge  of  the  oil  pocket, 
the  local  heat  flux  will  be  higher  than  the  average,  and 
hence  the  ice  will  grow  more  rapidly  there.  (The  heat  re- 
moved is  balanced  by  the  latent  heat  of  freezing  of  sea 

water.)  Before  the  oil  pocket  is  completely  included  in 
« 

sea  ice,  there  will  be  a rim  of  ice  around  the  pocket 
thicker  than  the  average.  It  may  be  that  the  existence  of 
such  a rim  might  be  an  aid  to  locating  oil  under  the  ice. 

(13) 

The  most  important  result  of  Wolfe  and  Hoult's 
work  is  an  estimate  of  how  much  oil  can  stick  to  the  under 
surface  of  the  sea  ice.  To  determine  this,  after  the  oil 
was  injected  into  the  freezing  apparatus,  the  ice  plug  was 
removed,  and  the  amount  of  oil  remaining  measured.  A care- 
ful study  was  made  of  the  various  mechanisms  by  which  a 
layer,  of' thickness,  6,  to  stick  to  the  under  surface  of 
the  ice. 

It  was  found  that  the  oil  adhered  to  the  ice  under  sur- 
face. If  Q/A  is  the  heat  flux  flowing  through  the  ice  and 
oil,  the  thickness,  6,  is  given  by  the  simple  formula f 


6-8  (Q/A) 

watts  ' 


In  the  Arctic  winter,  the  heat  flux  through  the  ice  is 
_2 

about  35  watts  m . This  gives  a 6 of  about  6.5mm.  If  crude 

oil  is  spilled  from  a supertanker,  with  a total  volume  of 
3 

113,000  m , the  total  area  which  could  be  covered  by  oil  is 
2 

about  17  km  (6  square  miles). 


This  is  an  upper  bound,  because  undoubtedly  some  of 
the  oil  will  be  trapped  In  deeper  pools. 


SUMMARY 

It  Is  helpful  to  conclude  this  report  by  commenting  on  the 
practical  application  of  the  results  of  this  work.  The  most  im- 
portant result,  which  is  substantiated  by  all  the  research  to 
date,  is  that  the  oil  in  even  a very  large  Arctic  spill  will  be 
confined  to  a very  small  area  due  to  natural  processes . This 
simple,  general,  conclusion  is  supported  by  both  the  spreading 
and  pocketing  experiments  as  well  as  the  studies  on  oil  evolu- 
tion. 

As  a corollary  to  this  result,  it  seems  clear  that,  unlike 
the  situation  in  temperate  waters,  there  is  little  incentive  to 
promptly  clean  up  large  oil  spills,  as  the  area  of  the  oil  spill 
will  not  measurably  increase  with  time,  while  the  oil  remains 
in  contact  with  sea  ice. 
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APPENDIX  A 


EFFECTS  OF  OIL  UNDER  SEA  ICE 


1.  INTRODUCTION 

Hie  discovery  of  oil  on  the  Northern  Slope  of  Alaska  had 
raised  many  questions  concerning  the  effects  of  oil  spills  In 
the  Arctic  Ocean  and  other  Ice  filled  waters.  The  handling 
and  transporting  of  oil  In  sizable  quantities  always  results 
in  some  spillage,  the  average  spillage  being  on  the  order  of 
0.12  of  the  quantity  transported  (Blumer,  1969,  p.  6).  The 
North  Slope  borders  on  the  Arctic  Ocean  which  is  essentially 
ice  locked  for  about  nine  months  of  the  year  (SEV  Data  Package , 
1970,  p.  11).  Even  during  the  summer  months  the  permanent  ice 
cap  remains  within  ten  nautical  miles  of  the  shore  (SEV  Data 
Package , 1970,  p.  34).  Any  sort  of  waterbom  transport  of 
North  Slope  crude  oil  will  create  the  possibility  of  spilling 
large  amounts  of  crude  oil  under  the  ice. 

Granted  that  under  ice  oil  spills  are  almost  certain  to 
occur  if  Arctic  oil  operations  continue,  it  is  Important  to 
know  how  such  spills  can  be  contained  and  removed  if  the  en- 
vironmental damage  caused  by  them  would  be  serious  enough  to 
warrant  such  action.  In  the  search  for  this  information  the 
first  question  to  ask  is  how  the  oil  spill  behaves  under  the 
ice.  Much  is  known  of  how  oil  spreads  on  temperate  waters, 
both  from  experimental  and  theoretical  studies  and  from  observa- 
tions of  large  scale  disasters  such  as  the  Torrey  Canyon  and 
Santa  Barbara  Channel  incidents.  Almost  nothing  is  known  of 
the  behavior  of  an  oil  spill  under  sea  ice.  The  speed  with 
which  the  oil  spreads  or  the  ultimate  thickness  to  which  it 
spreads  are  unknown.  The  interaction  between  the  oil  and  the 
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peculiar  microstructure  of  the  lower  layer  of  the  sea  ice  (see 
Weeks,  1966,  pp.  173-178)  is  also  unknown. 

There  has  been  a good  deal  of  speculation  concerning  the 
way  sea  ice  will  interact  with  trapped  oil  under  it.  Three 
principal  modes  of  behavior  have  been  considered  possible. 

(1)  The  sea  ice  will  entrap  the  oil  causing  the  formation  of 
a matrix  of  oil  and  ice . 

(2)  The  ice  will  entrap  the  oil  in  a pool  and  proceed  to  form 
beneath  it. 

(3)  The  ice  will  continue  to  grow  pushing  the  oil  before  it. 

(See  Fig.  1).  Such  a phenomenon  is  difficult  to  observe  in 
field  experiments,  so  laboratory  experiments  were  conducted  to 
determine  which  of  these  phenomena  occurred  and  to  examine  the 
phenomenon  quantitatively,  if  possible. 


EXPERIMENTAL  APPARATUS  AND  METHODS 

The  principal  function  of  the  laboratory  apparatus  was  to 
produce  a nearly  uniform  vertical  heat  flux  in  a tank  of  sea- 
water. The  apparatus  consisted  of  a 1.59  cm  (5/8  inch)  thick 
plexlglas  tank,  30. A cm  (12  inch)  square  in  cross  section  and 
1.06  m (3  1/2  feet)  deep  topped  by  a 30. A cm  square  stainless 
steel  cold  plate.  The  plate  was  connected  by  flexible  tubing 
to  a commercial  refrigerating  unit  which  circulated  freon  12 
refrigerant  at  a set  temperature  (+  1.1  deg)  dovm  to  -29°C. 

The  tank  was  tapered  0.03A  cm/cm  to  facilitate  removal  of  the 
ice  and  was  clad  on  all  sides  with  30. A cm  (12  inch)  thick 
urethane  foam  insulation.  The  insulation  was  installed  in 
separate  panels  which  could  be  removed  to  make  visual  observa- 
tions of  the  tank.  (See  Fig.  2).  As  a reasonable  facsimile 


FIGURE  2.  SEA  ICE  TEST  CELL  SCHEMATIC. 


of  sea  Ice  has  been  made  in  a cylinder  of  14  cm  diameter, 

(Lofgren  and  Weeks,  1969,  p.  157)  the  cross  section  of  the 
tank  was  not  a limiting  consideration.  The  convection  patterns 
were  presumed  sufficiently  small  that  proximity  of  the  wall  had 
no  effect  on  circulation  in  the  center  of  the  tank.  At  first 
an  attempt  was  made  to  make  the  heat  flux  through  the  insulating 
walls  very  small  in  comparison  with  the  heat  flux  through  the 
ice  when  producing  ice  at  Arctic  rates.  This  was  not  possible 
with  solid  insulating  materials  commercially  available,  so  to 
reduce  the  effects  of  heat  gain,  ice  was  produced  at  somewhat 
faster  rates  than  is  typical  for  Arctic  conditions  in  late 
winter. 

The  quantities  measured  during  the  experiments  were  the 
temperature  distribution  in  the  ice  and  water,  the  salinity 
of  the  water  beneath  the  ice,  the  depth  of  growth  of  the  ice 
and  the  amount  of  oil  injected  under  the  ice.  Temperature  was 
measured  with  a bridge  of  copper-constantan  thermocouples. 

The  salinity  of  the  water  beneath  the  ice  was  measured  by  the 
use  of  a specially  designed,  totally  submerged  hydrometer. 

The  depth  of  growth  of  the  ice  was  determined  by  measuring  the 
depth  as  seen  through  the  tank  with  a ruler.  The  growth  ve- 
locity of  the  ice  was  computed  by  noting  the  time  that  each 
thickness  measurement  was  made  and  calculating  the  arithmetic 
mean  velocity  between  successive  measurements. 

As  a substitute  for  seawater  solar  salt  was  dissolved  in 
tap  water  until  a salinity  of  approximately  30  o/oo  was  achieved. 
This  was  used  at  the  start  of  each  experiment , but  salt  rejec- 
tion by  the  ice  increased  the  salinity  of  the  water  to  between 
35  and  40  o/oo  by  the  time  oil  was  placed  under  it.  Two  types 
of  oil  were  used.  No.  2 diesel  was  chosen  because  it  is  readily 
available  and  because  it  is  a common  marine  fuel  which  may  be 
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subject  to  large  scale  spills.  The  other  oil  used  was  a North 
Slope  crude  oil. 

Each  experiment  was  started  with  an  initial  plate  temper- 
ature of  -29°C  until  a few  centimeters  of  ice  formed.  Then  the 

temperature  was  set  to  the  temperature  at  which  the  experiment 

was  to  run.  When  the  ice  reached  a depth  between  12  and  16  cm, 
oil  was  injected  through  the  bottom  of  the  tank  from  a plenum 
chamber  using  compressed  air  to  supply  the  driving  pressure. 
Except  in  a few  trial  runs,  enough  oil  was  added  to  completely 
cover  the  lower  surface  of  the  ice.  The  average  thickness  of 
the  oil  layer  ranged  between  1 and  2.6  cm.  The  oil  was  allowed 

to  remain  until  some  ice  began  to  form  on  its  lower  surface. 

This  took  12  - 24  hours  in  most  cases.  The  experiment  was  then 
terminated  and  the  ice  cube  removed. 

In  order  to  determine  the  extent  to  which  oil  became  en- 
trapped in  the  ice  or  adhered  to  its  lower  surface,  the  lower- 
most 2.5  cm  of  the  ice  block  was  sawed  off,  and  the  slab  of 
ice  was  melted.  The  dimensions  of  the  slab  were  noted  before 
melting.  The  oil  and  water  obtained  were  collected  in  a grad- 
uated cylinder  and  measured.  Because  the  technique  for  measure- 
ment of  the  adhering  oil  is  so  crude , the  limits  of  uncertainty 
for  this  data  are  rather  wide. 

Since  the  porosity  of  sea  ice  increases  markedly  at  tem- 
peratures near  the  melting  point,  an  attempt  was  made  to  see 
if  there  was  any  marked  change  in  the  mode  of  oil  entrapment 
as  the  ice  melts.  The  existing  apparatus  could  not  actually 
simulate  the  melting  of  sea  ice  as  it  occurs  in  the  ocean  be- 
cause the  tank  had  a finite  heat  gain  from  the  laboratory 
although  there  is  no  similar  source  of  heat  in  the  sea.  Never- 
theless, the  experiment  was  conducted.  After  the  completion 
of  an  experiment  performed  in  the  usual  manner,  the  temperature 


of  Che  cold  plate  was  raised  to  -1.4  C,  and  the  ice  was  allowed 
to  melt.  Temperature  profiles  were  recorded  every  two  hours 
and  visual  observations  of  the  Ice  were  made. 

A more  detailed  description  of  the  experimental  apparatus 
methods  Is  available  in  Wolfe,  1972,  pp.  14-22. 

3.  EXPERIMENTAL  OBSERVATIONS  AND  RESULTS 
3. 1 Sea  Ice  Growth  Rate 

A steady  state  approximation  to  the  flow  of  heat  from 
the  sea  water  to  the  cold  plate  may  be  used.  The  rate  of 
formation  of  ice  equal  to  the  rate  of  heat  flux  divided 
by  the  heat  of  formation  of  the  ice,  as  shown  below. 


Q/A 

piceL 


(1) 


Since  the  thermal  conductivity  of  the  ice  is  essentially 
constant  over  the  range  of  temperatures  encountered,  and 
the  flow  of  heat  is  uniform  upwards , the  temperature 

v 

gradient  in  the  ice  is  constant.  Since  the  cold  ice  sur- 
face is  at  the  top  of  the  sea  water,  the  "sea"  is  well 
mixed  and  of  uniform  temperature  approximately  equal  to 
the  liquidus  temperature  of  a NaCL  and  water  solution  of 
concentration  35  o/oo. 

The  temperature  data  can  be  nondimensionalized  and 
normalized  by  dividing  the  temperature  difference  between 
the  local  temperature  and  the  water  temperature  by  the 
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total  temperature  difference  between  the  water  and  the 
cold  plate. 


T - T 


T - 


plate 


- T 


Distance  is  normalized  by  measuring  the  distance  from  the 
lower  surface  of  the  ice  to  the  thermocouple  and  dividing 
this  by  the  total  thickness  of  the  ice. 


* 


z 


z/h 


All  of  the  temperature  data  obtained  before  the  addition 
of  the  oil  is  shown  plotted  in  this  manner  in  Fig.  3.  It 
demonstrates  that  within  the  limited  of  experimental  error, 
the  temperature  gradients  are  uniform  and  that  the  resistance 
of  the  thermal  boundary  layer  at  the  ice-water  interface  is 
negligible. 

The  heat  flow  from  the  warm  laboratory  to  the  test  cell 
retards  the  rate  of  formation  of  ice  below  that  predicted  by 
Eq.  (1)  above.  The  rate  of  heat  gain  by  the  ice  tank  from 
the  lab  should  be  nearly  constant  for  all  growth  conditions, 
and  thus  the  rate  of  formation  should  approach  that  given 
by  Eq.  (1)  for  high  rates  of  growth.  The  available  temper- 
ature and  growth  rate  data  may  be  nondimensionalized  and 
normalized  by  dividing  the  product  of  the  heat  of  forma- 
tion of  ice  and  the  measured  growth  rate,  pLV,  by  the  heat 
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FIGURE  3.  DIMENSIONLESS  TEMPERATURE  PROFILES. 


flux  out  through  the  Ice  as  predicted  by  Eq.  (2). 


Q/A  - 


k,  (AT). 
ice  Ice 

h. 

ice 


(2) 


This  dimensionless  growth  rate  is  plotted  as  a function  of 
the  heat  of  formation  and  growth  velocity  product,  pLV, 
divided  by  the  calculated  heat  gain  per  unit  area  through 
the  walls  of  the  tank  and  is  displayed  in  Fig.  4.  Note 
that  the  value  of  the  dimensionless  parameter  pLV/ [k(AT) /h] 
drops  noticeably  below  unity  when  the  ratio  pLV/ (0^a^n/A) 
is  about  150/1.  The  area  of  the  tank  walls  is  about  15 
times  the  area  of  the  cold  plate.  As  expected,  the  growth 
velocity  deviates  significantly  from  Eq.  (1)  when  the  ratio 
of  pLV  to  the  total  heat  gain  is  abouf  1/10.  The  large 
error  bars  on  each  data  point  are  due  primarily  to  the 
large  uncertainty  in  the  measurement  of  the  rate  of  growth. 

3.2  Mode  of  Oil  Entrapment 

The  actual  method  by  which  the  crude  and  diesel  oils 
are  trapped  in  the  sea  ice  is  a combination  or  modes  (A) 
and  (C)  shown  in  Fig.  1.  The  bulk  of  the  oil  is  pocketed 
in  a pool  below  the  original  ice  sub-surface  while  more 
sea  ice  proceeds  to  grow  under  it.  A small  amount  of 
oil  does  rise  into  the  pores  of  the  skeletal  layer,  and 
into  small,  vertical  shafts  which  rise  from  the  lower 
surface  of  the  ice.  These  shafts  may  be  air  bubbles  sim- 
ilar to  those  that  form  in  fresh  water  ice  due  to  the  in- 
ability of  air  entrained  in  the  water  and  rejected  upon 
freezing  to  escape  beyond  the  advancing  ice  front  (Weeks 


and  Assur,  1969,  p.  8).  The  total  amount  of  oil  contained 
In  these  air  bubbles  Is  small  In  comparison  to  the  amount 
of  oil  that  adheres  to  the  lowermost  inch  of  the  ice  (see 
below).  During  the  ice  melting  experiment,  a considerably 
greater  amount  of  oil  rose  into  the  ice,  especially  into 
the  air  shafts.  In  some  Instances  it  reached  as  high  as 
8 cm  above  the  lower  surface  of  a 13  cm  thick  block,  yet 
the  volume  of  oil  trapped  in  this  manner  was  still  small, 
amounting  to  no  more  than  1/150^  of  the  volume  of  the  ice. 

The  presence  of  the  oil  pool  has  a significant  effect 
upon  the  heat  transfer  through  the  ice.  The  simple  model 
of  heat  transfer  through  sea  ice  described  above  is  a uni- 
form solid  of  constant  conductivity  between  two  constant 
temperature  regions,  one  at  the  liquidus  temperature  of 
seawater  and  the  other  at  the  ambient  temperature  of  the 
air,  or  in  the  experimental  case,  the  cold  plate.  As  the 
oil  has  a thermal  conductivity  that  is  considerably  lower 
than  that  of  the  ice,  its  presence  acts  as  an  insulating 
layer,  impeding  the  flow  of  heat  and  reducing  the  temper- 
ature drop  across  the  ice. 

The  temperature  data  obtained  after  the  oil  was  placed 
under  the  ice  may  be  nondimenslonalized  in  the  same  manner 
as  was  done  for  the  data  taken  with  ice  only.  The  origin 
for  the  vertical  axis  remains  at  the  lower  surface  of  the 
original  ice  block.  As  a result  temperature  measurements 
in  the  oil  and  below  are  expressed  as  negative  numbers. 

The  thickness  of  the  ice  before  the  oil  was  added  is  used 
to  normalize  distance,  but  the  total  temperature  drop 
across  the  ice  and  oil  is  used  to  normalize  temperature. 

The  normalized  temperature  distribution  are  shown  plotted 
in  Fig.  5 for  all  temperature  measurements  taken  after 
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FIGURE  5.  DIMENSIONLESS  TEMPERATURE  DISTRIBUTION  AFTER 
OIL  INSERTION. 


the  oil  was  in  place. 

As  can  be  seen  in  Fig . 5 , the  temperature  in  the  ice 
is  linear  at  all  times  after  the  addition  of  the  oil  ex- 
cept for  the  case  of  Experiment  9,  in  which  the  tempera- 
ture of  the  plate  was  high  enough  that  the  heat  gain  by 
the  tank  caused  the  water  temperature  to  rise  continuously. 
In  the  four  other  cases  in  which  temperature  profiles 
were  measured,  the  temperature  gradient  in  the  ice  remained 
constant,  but  decreased  in  magnitude  from  its  initial  value 
as  shown  in  Fig.  3 through  intermediate  values  to  the  final 
value.  The  temperature  distribution  in  the  oil  could  not 
be  determined  directly,  since  the  oil  layer  was  between 
1.3  and  2.5  cm  thick  and  the  thermocouples  were  spaced 
2.54  cm  (one  inch)  apart  in  the  bridge.  It  is  clear  that 
temperature  gradients  in  the  oil  were  much  larger  than 
in  the  ice. 

Whereas  clearly  the  heat  transfer  through  the  ice  is 
that  of  solid  body  conduction,  the  mode  of  heat  transfer 
in  the  oil  pool  is  not  obvious.  It  should  be  noted  in 
Fig.  5 that  for  all  of  the  crude  oil  experiments,  the 
ratio  of  the  temperature  drop  across  the  oil  to  that 
across  the  oil  and  ice  is  the  same,  0.6.  This  ratio  is 
independent  of  the  thickness  of  the  oil  pool,  which  varies 
between  1.3  and  2.5  cm  (in  dimensionless  units,  (z/h)  , 
-0.097  and  -0.189).  For  the  diesel  oil  the  ratio  of  the 
temperature  drop  across  the  oil  to  the  total  temperature 
drop  across  the  oil  and  ice  is  0.3.  This  behavior  is  not 
in  accordance  with  the  theory  of  static  heat  conduction. 

The  thermal  conductivity  of  most  petroleum  distillates  is 
only  a very  weak  function  of  temperature,  and  over  the 
range  of  temperatures  considered,  the  thermal  conductivity 
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of  the  liquid  fractions  of  petroleum  changes  less  than  3%. 
(Wolfe,  1972,  Appendix  I).  Within  the  limits  of  precision 
of  tl.  *.s  experiment,  it  may  be  considered  constant.  For 
constant  thermal  conductivity,  the  thermal  resistance  of 
a solid  conducting  heat  in  one  direction  is  h/k,  and  thus 
the  temperature  drop  across  the  solid  is  proportional  to 
its  thickness.  For  two  solids  conducting  in  series,  the 
ratio  of  the  temperature  drop  across  one  to  the  temperature 
drop  across  the  other  is  ^2^^)  • This  is  clearly 

not  the  case  for  the  ice  block  and  oil  pool. 

It  is  not  proper  to  assume  that  the  mode  of  heat 
transfer  in  the  oil  pool  is  by  conduction  alone.  At  the 
temperature  involved,  the  crude  oil  is  indeed  highly 
viscous,  but  further  examination  is  necessary  to  deter- 
mine if  it  may  be  considered  a solid  body.  The  diesel 
fuel  is  quite  fluid;  it  is  entirely  possible  that  con- 
vection cells  could  form  and  reduce  the  thermal  resistance 
of  the  oil  layer  below  that  predicted  by  static  conduction 
theory.  The  Nusselt  number  is  the  dimensionless  ratio 
which  compares  the  ratio  of  the  convective  heat  transfer 
coefficient,  K,  to  the  conductive  heat  transfer  coefficient, 
k/h.  For  free  convection  between  plane,  horizontal  sur- 
faces, the  Nusselt  modulus  may  be  defined  as 


Nu  = K h/k 


(3) 


The  Nusselt  modulus  can  be  shown  to  be  a function 
of  the  Grashof  number  and  the  Prandtl  number,  defined 
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by  relations, 


Or  = a^f.V 


(A) 


Pr  = c y/k 


(5) 


provided  that  the  mode  of  heat  transfer  is  free  convec- 
tion only  (Kreith,  1965,  p.  333).  When  a normal  fluid 
whose  density  decreases  with  increasing  temperature  is 
placed  between  plane  surfaces  and  cooled  from  above, 
convection  results.  This  convection  is  commonly  described 
in  terms  of  the  Rayleigh  number  which  is  the  product  of 
the  Grashof  and  Prandtl  numbers,  GrPr.  When  this  dimen- 
sionless variable  is  less  than  about  1700,  no  convection  re 
suits,  and  heat  is  transferred  by  static  conduction  only. 
When  GrPr  is  between  1700  and  42,000  convection  is  in  the 
peculiar  cellular  manner  known  as  Benard  cell  convection. 
Irregular  turbulence  results  at  values  of  GrPr  Higher  than 
42,000  (Eckert  and  Drake,  1959,  p.  328).  For  values  of 
GrPr  less  than  1700,  the  Nusselt  modulus  is  by  definition, 
unity.  Benard  cell  convection  may  be  adequately  correlated 
as  a function  of  Nusselt  modulus  by  Eq.  (6) , (Eckert  and 
Drake,  1959,  p.  328). 


Nu  - 0.107  (GrPr) 


0.3 


(6) 
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The  Nusselt  number  for  the  oil  pool,  can  be  expressed 
In  terras  of  the  measured  parameters  of  the  experiment. 

For  one  dimensional  heat  flow,  the  ratio  of  the  thermal 
resistances  of  the  oil  to  the  ice  is  equal  to  the  ratio 
of  the  temperature  drops  across  the  ice  and  across  the 
oil,  respectively.  Since  the  value  of  the  thermal  re- 
sistance of  the  oil,  is  expressible  in  terms  of 

known  parameters,  (AT)icg,  kicg , the  Nusselt  num- 

ber of  Eq.  (3)  may  be  evaluated  by  Eq.  (7). 


Nu  - 


k.  h (AT). 
ice  oil  7 ice 

k ..  h.  (AT)  .. 
oil  ice  ' 7 oil 


The  Rayleigh  number  shown  above  has  been  evaluated  for 
each  of  the  experiments  in  which  temperature  profiles  were 
measured.  The  Rayleigh  number  for  the  diesel  oil  measure- 
ment was  about  22,000  while  GrPr  for  the  crude  oil  data 
ranged  between  770  and  9700.  A log-log  plot  of  the  Nusselt 
number  as  a function  of  Rayleigh  number  is  shown  in  Fig.  6. 
The  data  conforms  reasonably  well  to  the  standard  relation- 
ship of  Eq.  (6)  which  is  sketched.  Although  the  data  could 
be  fitted  better  with  another  power  law,  it  would  be  unwise 
to  do  so,  since  so  few  measurements  have  been  made.  The 
deviation  from  Eq.  (6)  may  simply  be  due  to  errors  in 
measurement  since  the  data  from  which  it  was  derived  ex- 
hibits substantial  scatter  for  relatively  low  Rayleigh 
numbers  (Eckert  and  Drake,  1959,  p.  331).  Furthermore, 
Nusselt  numbers  less  than  unity  must  be  attributable  to 
errors  in  measurement. 

As  temperature  measurements  could  not  be  made  for 
all  experiments,  the  heat  flux  through  the  ice  after  the 
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FIGURE  6.  NUSSELT  MODULUS  AS  A FUNCTION  OF  RAYLEIGH 
NUMBER  FOR  BENARD  CELL  CONVECTION  IN  THE 
OIL  POOL. 


oil  Is  in  place  must  be  estimated  from  the  temperature  of 
the  cold  plate  only.  It  was  observed  that  the  temperature 
drops  across  the  oil  and  ice  were  constant  fractions  of 
the  total  temperature  drop  between  the  cold  plate  and  the 
sea  water,  irrespective  of  the  thickness  of  the  ice  or  of 
the  oil  but  dependent  upon  the  type  of  oil.  That  this  is 
expected  can  be  shown  from  the  equations  for  one  dimen- 
sional heat  conduction.  The  expressions  for  the  flow  of 
heat  through  the  ice  and  oil  can  be  written  in  terms  of 
the  thermal  resistences  and  temperature  drops  in  the 
system. 


(AT) 


ice 


(AT) 


ice 

lice 


total 

1 

' + K 


(8) 


The  linear  temperature  gradient  in  the  ice  is  expressible 
in  terms  of  the  Nusselt  number  by  substituting  Eq.  (3) 
into  Eq.  (3)  and  rearranging. 


(AT) 


(AT) 


total 


ice 


1 + 


kjcc  ^oil  -1 

k h.  Nu 
oil  ice 


The  expression  for  the  Rayleigh  number  is 


GrPr 


cp3  g3  (AT)  ho3 


(9) 


(10) 
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For  constant  properties  of  the  oil  and  for  small  variations 
in  (AT),  the  Rayleigh  number  is  pronortional  to  the  cube  of 
the  oil  thickness.  Equation  (6)  shows  the  Nusselt  number 
is  proportional  to  the  0.3  power  of  OrPr.  Therefore,  for 
small  variations  in  (AT)  and  in  the  fluid  properties,  the 
approximate  equality  Nu  - Constant  , (hQ^)  is  valid. 
Substituting  this  in  Eq.  (10) 


(AT) 


ice 


1 + 


kice 


<koil> 


ice 


Constant 


(ID 


It  can  be  seen  that  for  small  variations  in  h.  , (AT). 

ice’  ' 'ice 

is  a nearly  constant  fraction  of  The  thickness 

of  the  ice  in  the  experiments  ranged  between  13  and  17  cm. 

It  is  therefore  reasonable  to  consider  that  the  ratio  of 
the  temperature  drops  (AT)^ce/(AT)  is  the  same  for  all 
experiments  as  it  was  for  those  in  which  the  temperature 
profiles  were  measured.  In  all  subsequent  calculations 
involving  the  heat  flux  through  the  ice  after  the  oil  is 
present,  (AT)icg/(AT)  will  be  taken  to  be  0.4  for  the 
crude  oil  experiments  and  0.7  for  the  diesel  oil  experiments. 


3 . 3 Thickness  of  Oil  Adhering  to  the  Ice 

Upon  removal  of  the  ice  block  from  the  tank,  some  of 
the  oil  remained  with  the  ice,  either  trapped  in  the  porous 
ice  structure  or  adhering  to  the  lower  surface,  while  the 
bulk  of  the  oil  remained  in  the  original  oil  pool.  The 
coating  of  diesel  oil  was  fairly  uniform  but,  the  crude 
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oil  was  quite  uneven  In  Its  adherence  to  the  Ice.  It  some- 
times coated  the  Ice  In  a continuous  layer  of  uneven  thick- 
ness, but  in  other  cases  it  separated  from  the  ice  in 
patches,  which  were  clean  of  oil  (except  for  the  oil  in 
the  interstices  of  the  sketetal  layer).  A photograph  of 
an  ice  block  after  removal  from  the  tank  following  a 
typical  crude  oil  experiment  is  shown  in  Fig.  7.  Although 
the  oil  does  not  adhere  to  the  ice  in  an  even  layer,  an 
equivalent  thickness  of  the  oil  layer  may  be  calculated 
by  dividing  the  total  volume  of  oil  measured  after  melt- 
ing the  ice  sample  by  the  area  of  the  subsurface  of  the  ice 
sample.  It  wa3  difficult  to  determine  visually  to  what 
extent  the  oil  was  traooed  in  the  interstices  of  the  ice 
and  to  what  extent  it  simply  adhered  to  the  outside.  Most 
of  the  oil  appeared  to  be  confined  to  the  lowermost  0.5  cm, 
although  a section  of  ice  this  thin  could  not  be  sliced, 
it  is  certain  that  all  but  slight  traces  of  the  oil  was 
contained  within  the  lowermost  2.5  cm  of  the  ice.  The 
traces  above  this  level  were  due  exclusively  to  the  presence 
of  the  air  bubbles . 

The  details  of  the  mechanism  by  which  oil  adheres  to 
the  ice  are  neither  simple  nor  obvious.  From  physical 
considerations  it  is  possible  to  infer  that  the  process 
of  oil  adhesion  is  both  local  and  steady  state.  The  ad- 
hesion phenomenon  can  not  reasonably  be  expected  to  de- 
pend directly  on  any  properties  of  the  ice  far  away  from 
the  ice-oil  interface.  Thus  the  thickness  of  the  ice  or 
of  the  oil  pool,  the  absolute  temperature  of  the  cold 
plate,  the  rate  of  growth  of  ice  below  the  oil  pool  or  the 
flow  field  within  the  oil  pool  can  not  be  expected  to  in- 
fluence the  thickness  of  the  adhering  layer  of  oil,  6. 
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FIGURE  7.  PHOTOGRAPH  OF  AN  ICE  BLOCK 
UPON  REMOVAL  FROM  THE  TEST 
CELL  AFTER  A TYPICAL  CRUDE 
OIL  EXPERIMENT. 
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The  oil  thickness  may  depend  upon  the  geometry  of  the  ice 
subsurface,  the  bulk  properties  of  the  ice  and  oil  at  the 
interface,  the  surface  energy  of  the  interface,  and  the 
heat  flux  through  the  interface. 

That  the  adhesion  process  is  steady  state  may  be  in- 
ferred from  the  fact  that  there  were  no  transient  effects 
taking  place  upon  removal  of  the  ice  from  the  tank.  At 
this  time,  all  temperatures  were  in  accordance  with  those 
predicted  by  static,  one  dimensional  heat  conduction. 

There  was  no  discernable  advance  of  the  ice  front  after  the 
oil  was  placed  under  the  ice,  and  there  was  no  observable 
change  in  the  amount  of  oil  entrapped  after  the  steady 
state  temperatures  were  reached.  There  was  no  trend  in  the 
adhesion  thickness  measured  as  a function  of  the  time  it 
remained  under  the  ice. 

The  pronerties  of  the  ice  and  oil  at  the  interface 
which  might  be  of  significance  are  the  viscosity  of  the 
oil,  the  buoyancy  of  the  oil,  the  mechanical  properties  of 
the  ice  and  the  interfacial  energy  of  the  ice  and  oil. 

Since  the  ice  does  not  fracture  or  deform  as  a result  of  the 
oil  adhesion,  its  mechanical  properties  could  be  of  no  ef- 
fect on  6.  Since  the  oil  is  not  in  motion  when  the  lower 
surface  is  removed,  the  viscosity  of  the  oil  could  not  be 
of  influence  on  6.  The  heat  flux  through  the  interface  is 
a function  of  the  steady  state  thermal  conductivities  of 
the  ice  and  oil,  but  not  of  the  thermal  diffusivity  which 
specifies  the  transient  conduction  properties  of  the  ma- 
terials. The  interfacial  energies  of  the  ice  and  oil  can 
be  described  in  terms  of  a bulk  property  known  as  the  ad- 
hesion coefficient. 
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In  the  simple  theory  of  adhesion,  the  work  necessary 
to  create  an  additional  unit  surface  may  be  expressed  in 
terms  of  an  adhesion  coefficient,  w,  which  has  the  same 
dimensions  as  surface  tension.  The  simplest  theory  for 
adhesion  coefficients  relates  them  directly  to  surface 
tension,  0. 


u ■ aa/v  + - °./s 


(12) 


in  which  a , and  a , are  the  common  surface  tensions 
a/v  s/v 

usually  measured  for  liquids  or  solids,  i.e.  the  energy 

needed  to  create  a unit  surface  in  a vacuum,  and  the 

term  a , is  a measure  of  the  interfacial  energy  per  unit 
a/s 

surface  of  interface  between  the  adhesive  and  the  substrate 
(Skeist",  1962,  p.  38).  Although  such  quantities  are  simple 
to  state  in  principle,  they  are  not  easy  to  measure.  How- 
ever, it  is  true  for  both  water  and  petroleum  derivatives 
that  the  surface  tensions  do  not  vary  appreciably  over  the 
range  of  temperatures  encountered  in  this  experiment. 
Therefore,  a preliminary  estimate  of  the  surface  tensions 
may  be  found  by  using  the  tabulated  values  for  the  surface 
tension  between  oil  and  air,  and  water  and  air  for 
and  respectively,  and  the  tabulated  values  for  the 

interfacial  tension  between  oil  and  water  for  Ofl^g. 

(Wolfe,  1972,  Appendix  I).  Whether  this  value  is  accurate 
or  not  does  not  preclude  its  use  in  evaluating  the  im- 
portance of  the  adhesion  coefficient.  Surface  tension  is 
a strong  function  of  the  materials  interfacing  each  other 
and  only  a weak  function  of  the  local  temperature.  Care 
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was  taken  to  insure  that  the  sea  Ice  and  oil  used  in  the 
experiments  were  of  the  same  composition  as  they  would 
be  in  the  Arctic.  Thus,  no  matter  what  the  actual  value 
of  w,  it  can  be  expected  to  be  the  same  value  in  the  field 
that  it  was  in  the  laboratory. 

Elementary  physical  considerations  have  reduced  the 
list  of  parameters  which  may  be  of  influence  on  5 to  the 
buoyancy,  pAg,  the  adhesion  coefficient,  w,  the  heat  flux 
through  the  ice,  0/A,  and  the  geometric  properties  of  the 
Interface  characterized  by  the  intersticial  spacing  of 
the  dendrites.  (Weeks,  1966,  p.  175).  The  intersticial 
spacing  is  a function  of  the  velocity  of  growth  of  the  ice 
before  the  oil  is  in  place  and  is  specified  by  Eq.  (13) 


aQ  - 3.95  x 10“5/»^dh7dt  (SI  units)  (13) 

(Assur  and  Weeks,  1963,  pp.  97-98).  Since  the  diameter  of 
the  pores  of  the  skeletal  layer  is  roughly  constant,  the 
volume  of  the  skeletal  layer  filled  with  water  can  be  ex- 
pected to  get  larger  as  aQ  becomes  smaller.  Any  entrap- 
ment phenomenon  that  is  dependent  on  filling  the  pores  of 
the  skeletal  layer  with  oil  would  be  expected  to  produce 
larger  values  of  6,  as  the  growth  velocity  of  the  ice  be- 
fore the  oil  is  added  becomes  larger.  In  fact  just  the 
opposite  occurred.  Less  oil  adhered  at  higher  growth 
velocities . 

Dimensional  analysis  can  be  used  to  determine  the 
relationship  between  the  heat  flux  through  the  ice,  the 
buoyant  forces  on  the  oil  and  the  adhesive  forces  on  the 
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oil.  If  the  thickness  of  the  adhering  oil  is  presumed  to 

be  a function  of  heat  flux  and  buoyancy  only,  the  Buckingham 

Pi  Theorem  indicates  that  only  one  dimensionless  grouping 

can  be  formed,  and  that  is  invariant  for  all  values  of 

density  and  heat  flux.  That  dimensionless  variable  is 
2 2 1/3 

6[(pAg)  g/(Q/A)  ] . Since  density  is  the  same  in  the 

Arctic  as  in  the  laboratory,  the  experimental  values  of  6 
should  be  extrapolated  to  field  conditions  by  considering  6 
to  be  a linear  function  of  (Q/A)  . By  similar  reason- 

ing if  the  adhesion  thickness  is  considered  to  be  a func- 
tion of  the  adhesion  coefficient,  the  heat  flux  through 
the  ice  and  oil,  and  the  acceleration  of  gravity,  the  only 
dimensionless  grouping  possible  is  6 (Q/A)  /w  g,  and  it 
must  be  invariant.  Since  w is  the  same  in  the  Arctic  as 
in  the  experiment,  the  experimental  values  of  6 are  ex- 
pected  to  be  directly  proportional  to  (Q/A)  “. 

A dimensional  log-log  plot  of  the  adhesion  thickness, 

6,  as  a function  of  (Q/A)  was  made  in  order  to  reveal  which 
power  law  predominates.  (See  Fig,  8),  As  can  be  seen  from 
the  plot,  within  the  limits  of  orecision  of  the  experiment, 
the  data  conforms  to  a power  law  of  the  form 


6 - C(Q/A)"2 


(14) 


This  indicates  that  within  the  limits  of  precision  of  this 
experiment,  the  data  correlates  as  a function  of  the  di- 
mensionless  parameter  6 (Q/A)  /w  g,  and  is  therefore  de- 
pendent upon  the  heat  flux  through  the  ice  and  oil,  and 
the  value  of  the  adhesion  coefficient  alone.  The  limits 
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of  precision  of  the  experiment  are  not  very  good.  The 
coefficient,  C,  of  Eq.  (14)  may  be  evaluated  within  a 
factor  of  2 or  3 at  best.  The  "best"  power  law  approxi- 
mation to  the  data  shown  in  Fig.  8 was  drawn  freehand  and 

3-1 

yields  a value  of  C which  is  8 in  v . The  equations  for 

the  limits  of  uncertainty  of  the  data  are  also  of  the  form 

of  Eq.  (14)  and  yield  maximum  and  minimum  values  of  C which 
3 -1 

are  16  and  3 m w , respectively.  Equation  (14)  may  be  re- 
written in  dimensionless  form 

A - (14a) 

w g 

The  dimensionless  constant  A has  a best  value  of  330  and 
may  vary  between  110  and  660  due  to  the  uncertainty  of  the 
data.  The  adhesion  coefficient,  w,  has  been  taken  to  be 
0.05  n m 

Equation  (14)  may  be  used  to  extrapolate  the  measure- 
ments of  this  experiment  to  actual  Arctic  conditions  with 
an  error  factor  of  between  2 and  3.  As  shown  below,  such 
extrapolation  can  be  used  to  give  an  upper  bound  to  the 
extent  to  which  an  oil  spill  can  spread  in  the  Arctic. 

The  correlation  of  the  experimental  data  shown  above  does 
not  yield  a great  deal  of  insight  into  the  physical  mech- 
anism by  which  the  oil  adheres  to  the  ice.  It  does  indi- 
cate that  the  mechanism  is  related  to  the  interfacial 
energies  of  the  oil  and  the  ice,  rather  than  to  the  geometry 
of  the  surfaces  or  the  other  bulk  properties  of  the  oil. 

In  general,  the  phenomenon  of  adhesion  is  poorly  understood, 
and  a detailed  examination  of  the  actual  mechanism  involved 
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is  beyond  the  scope  of  this  work. 

It  remains  to  explain  the  fact  that  almost  no  oil 
was  observed  to  adhere  to  the  ice  during  one  of  the  crude 
oil  experiments . This  occurred  during  the  coldest  experi- 
ment run  with  a crude  oil  sample  in  which  the  steady  state 
temperature  gradient  in  the  ice  after  oil  injection  was 
86.6  deg  m and  the  mean  temperature  of  the  crude  oil  was 
-14.5°C.,  (well  below  the  pour  point  for  North  Slope  crude 
oil  of  approximately  -4°C) . Upon  removal  of  the  ice  block, 
the  oil  sheared  neatly  awav  from  the  ice  at  the  ice-oil 
interface.  The  oil  appeared  to  be  frozen,  and  was  the  con- 
sistency of  ice  cream.  It  may  be  that  when  the  oil  freezes, 
the  mechanism  of  oil  adhesion  is  altered  because  the  oil 
behaves  as  a solid,  instead  of  a liquid.  If  the  block  of 
oil  is  sufficiently  heavy,  its  weight  could  cause  it  to 
break  away  from  the  icc  at  the  ice-oil  interface,  because 
this  is  the  weakest  point  in  the  composite  ice  and  oil 
structure.  Such  an  observation  is  not  of  much  practical 
significance  since  temperature  gradients  of  this  magnitude 
rarely  occur  in  natural  sea  ice. 


4.  CONCLUSIONS  AND  PRACTICAL  APPLICATIONS 

Three  important  conclusions  can  be  drawn  from  these  experi- 
ments concerning  the  behavior  of  crude  oil  when  trapped  under 
sea  ice.  First,  the  mode  of  entrapment  of  the  oil  has  been 
conclusively  determined.  Second,  the  effect  of  the  presence  of 
an  oil  pool  upon  the  growth  rate  of  the  ice  and  upon  its  tem- 
perature distribution  has  been  determined.  Third,  an  order  of 
magnitude  estimate  of  the  maximum  extent  to  which  oil  can 
spread  under  Arctic  ice  can  be  deduced  from  the  measurements 
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of  the  thickness  of  the  adhering  oil  layer. 

The  experiments  have  conclusively  determined  the  mode  by 
which  oil  is  entrapped  under  a growing  sheet  of  ice.  The  ex- 
tent of  the  oil's  entrapment  in  the  ice-brine  matrix  has  been 
shown  to  be  negligible.  Even  when  the  porosity  of  the  ice  in- 
creases markedly  during  prcmelt  conditions,  the  volume  of  oil 
entrapped  remains  small.  Also,  it  has  been  shown  that  the  ice 
does  not  in  any  way  grow  through  the  oil.  The  oil  is  neatly 
pocketed  by  the  ice  as  more  ice  proceeds  to  form  under  it.  Tn 
the  absence  of  currents  under  the  ice,  even  the  large  pockets 
of  oil  which  might  form  in  the  larger  recesses  of  the  ice  sub- 
surface can  be  expected  to  be  entrapped  as  a whole.  Knowledge 
of  this  behavior  should  have  significant  impact  on  the  design 
of  equipment  and  procedures  to  clean  up  oil  spills  under  ice. 

The  presence  of  an  oil  pool  under  the  ice  causes  a marked 
change  in  the  temperature  distribution  in  the  ice  because  the 
oil  pool  acts  as  an  insulating  layer  between  the  cold  air  and 
the  relatively  warmer  sea  water.  The  temperature  distribution 
in  the  ice  can  be  calculated  by  a one  dimensional  analysis  for 
oil  pools  whose  widths  are  large  in  comparison  with  their  depths. 
The  exact  procedure  for  such  a calculation  would  depend  upon  the 
quantities  which  are  easily  measurable  in  a given  oil  spill. 

If,  for  example,  the  temperature  distribution  in  the  ice  above 
the  oil  were  desired,  and  the  known  parameters  were  the  thick- 
ness of  the  ice  before  the  spill,  h^ce^,  the  thickness  of  the 
oil  pool,  the  thickness  of  the  ice  under  the  oil  pool,  ^^ce2» 
and  the  temperature  of  the  upper  ice  surface,  then  by  presuming 
the  heat  flux  uniform,  the  temperature  difference  across  h^ce2 
could  be  calculated  by  the  following  relation 


(AT) 


(h,  ./ k.  ) (AT)  _ . 

icel  ice  total 


icel 


icel 

Cice 


1 lce2 

K k. 

— ice 


(15) 
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The  difficulty  with  this  procedure  is  that  the  thermal  resistance 
of  the  oil,  1/K,  is  dependent  upon  the  temperature  drop  across 
it  in  accordance  with  Eqs.  (3)  through  (6),  and  is  thus  dependent 
on  (AT)icg  and  (AT)  Since  Eq.  (15)  is  non-linear  and  not 

explicit  in  (AT)^cg  it  must  be  solved  by  some  sort  of  iterative 
procedure,  such  as  simple  trial  and  error,  but  it  can  be  solved, 
since  all  the  necessary  governing  relations  have  been  defined. 

The  distribution  of  temperature  in  the  ice  block  is  linear,  so 
the  temperature  gradient  may  be  readily  calculated  to  determine 
the  rate  of  growth  under  the  ice.  Also,  this  approach  may  be 
extended  to  cover  cases  in  which  additional  thermal  resistances, 
such  as  those  due  to  snow  cover  or  the  thermal  boundary  layer 
of  the  air,  are  present.  The  effects  of  radiative  heat  trans- 
fer from  the  ice  surface  also  may  be  included.  In  a converse 
procedure,  the  thickness  of  an  oil  pool  under  the  ice  may  be 
inferred  from  measuring  the  thickness  of  the  ice  and  the  tem- 
perature gradient  within  it.  Whether  such  calculations  are  of 
practical  import  remains  to  be  seen. 


The  experiments  which  measured  the  effective  thickness  of 
adhering  oil  layer,  6,  con  be  used  to  give  an  upper  bound  esti- 
mate of  the  extent  to  which  an  oi7  spill  under  ice  can  spread. 

It  can  be  shown  that  the  acceleration  applied  to  the  oil  layer 
by  currents  under  the  ice  is  much  smaller  than  the  acceleration 
due  to  gravity  applied  to  the  oil  layer  when  it  is  lifted  from 
the  tank  prior  to  slicing.  Hie  shear  stress  applied  to  a plane 
surface  in  turbulent  flow  can  be  expressed  in  terms  of  the  fluid 
density,  the  free  stream  velocity,  and  a friction  factor  which 
is  a function  of  the  surface  roughness  and  nearly  independent 
of  Reynolds  number.  (Rohsenow  and  Choi,  1961,  pp.  76-78). 


T - 


(16) 
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The  acceleration  applied  to  the  fluid  can  be  defined  in  terms 
of  the  friction  velocity,  U , as  follows 


8drag 


(17) 


Where  the  mean  roughness  height,  e,  is  the  approximate  eddy  size 
at  the  surface.  For  a roughness  height  of  0.5  m,  a current  ve- 
locity of  0.5  ra  s ^ (approximately  1 knot)  and  a friction  factor 

of  0.04,  the  effective  acceleration  is  calculated  by  Eqs.  (16) 

_2 

and  (17)  to  be  0.01  m s , which  is  very  small  in  comparison  with 
the  acceleration  due  to  gravity.  Thus,  attempts  to  calculate 
the  maximum  area  over  which  the  oil  spreads  by  using  the  effec- 
tive oil  thickness  measured  in  these  experiments  will  be  quite 
conservative . 

The  effective  thickness  of  the  oil  adhering  to  the  ice,  6, 
determines  the  extent  to  which  a given  volume  of  oil  will  spread 
according  to  the  relation 


A 


Volume  of  oil 

S 


(18) 


It  has  been  shown  that  6 is  a function  of  the  heat  flux  through 
the  ice  in  accordance  with  Eq.  (14).  Thus  from  knowledge  of  the 
growth  conditions  of  the  ice  and  the  volume  of  oil  spilled,  the 
area  over  which  the  oil  spreads  can  be  determined.  For  example, 
typical  mid-winter  conditions  for  the  Northern  Slope  would  be 
characterized  by  the  presence  of  first  year  sea  ice  approxi- 
mately 1.5  m thick  with  an  air  temperature  of  -15°c. 
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(SEV  Data  Package,  1970,  p.  11).  Ignoring  the  effects  of  con- 
vective and  radiative  heat  transfer  at  the  ice-air  boundary  and 

assuming  that  no  snow  cover  is  present,  the  heat  flux  through 

-2 

the  ice  can  be  calculated  to  be  approximately  35  watts  m . 

Since  this  calculation  is  performed  to  determine  the  maximum 
spread  of  the  oil,  the  thickness  of  the  oil,  and  thus  its  ther- 
mal resistance,  will  be  very  small  in  comparison  with  that  of 
the  1.5  m thick  ice.  The  heat  flux  through  the  oil  and  the  ice 

can  be  considered  to  be  the  same  as  for  the  ice  alone,  35  watts 
_2 

m . From  Eq.  (14),  6 is  calculated  to  be  6.5  mm  with  a maxi- 
mum value  of  13  mm  and  a minimum  of  2.5  mm.  If  crude  oil  spil- 
led is  from  a supertanker  of  100,000  metric  tons  capacity,  the 
total  volume  of  oil  will  be  113,000  m , and  from  Eq.  (18)  the 

2 

total  area  over  which  it  will  spread  will  be  a maximum  of  45  km 

2 

(17  square  miles),  a minimum  of  9 km  (4  square  miles)  and  a 

2 

"best"  estimate  of  17  km  (6  square  miles).  Even  the  largest 
of  these  values  is  considerably  smaller  than  the  area  over  which 
such  a spill  would  spread  in  temperate  waters.  Similar  calcula- 
tions can  be  performed  for  the  conditions  of  a particular  spill, 
in  accordance  with  the  local  ice  thickness  and  temperature  con- 
ditions . 

It  should  be  pointed  out  that  all  of  the  above  conclusions 
apply  to  conditions  typical  for  first  year  sea  ice.  Multi 
year  ice  is  generally  of  much  lower  salinity  and  different  com- 
position due  to  the  fact  that  brine  drains  from  it.  Because  it 
is  older,  its  history  is  less  certain  and  because  it  occurs 
primarily  in  the  permanent  ice  pack,  it  has  often  been  deformed 
and  distorted  by  the  large  forces  within  the  ice  cap.  Yet  from 
the  physical  analysis  of  the  oil  adhesion  phenomenon,  none  of 
the  above  differences  between  new  and  old  ice  should  affect  the 
ultimate  thickness  to  which  the  oil  can  spread.  In  mid-winter 


mam mmm 


the  lowermost  edge  of  the  old  Ice  will  be  composed  of  new  Ice 
which  hes  just  formed.  The  heat  flux  will  be  substantially 
lower  because  the  old  ice  is  thicker,  and  thus  Eq.  (14)  will 
predict  larger  values  of  6 than  would  occur  for  new  ice  under 
similar  weather  conditions. 
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APPENDIX  B 


EXPERIMENTAL  DETERMINATION  OF  THE  EXTENT  OF  AN 
OIL  SPILL  ON  ARCTIC  ICE 


1.  INTRODUCTION 

In  1968,  oil  was  discovered  on  the  North  Slope  of  Alaska. 
Since  the  world  supplies  of  oil  are  dwindling  and  the  prices 
rising,  this  oil  must  be  used  if  the  present  standard  of  living 
is  to  be  maintained  at  a reasonable  cost.  This  means  that  the 
crude  oil  must  be  shipped  to  a place  where  it  can  be  refined, 
and  whether  supertanker  or  pipeline  is  used,  the  oil  will  have 
to  travel  over  the  rough  surfac  s of  the  Arctic.  A small  frac- 
tion of  all  oil  shipped  is  spilled.  In  order  to  deal  with  this 
environmental  problem  it  appears  necessary  to  gain  an  understand- 
ing of  how  oil  spreads  on  or  under  ice.  This  understanding  must 
come  in  two  parts:  first,  the  process  of  spreading  must  be 

understood  and  the  time  for  spreading  determined;  and  second, 
the  final  size  of  the  oil  pool  must  be  determined. 

Initially,  an  inertial  spreading  law  was  proposed  to  explain 
the  first  problem.  Using  this  theory  the  rate  of  spread  could 
be  scaled  using  the  volume  released  and  a gravity  terra.  The 
spread  would  be  Independent  of  the  surface  characteristics.  In 
order  to  solve  the  second  problem,  one  must  take  into  account  the 
characteristics  of  the  ice  surface,  namely  its  roughness  and  the 
random  nature  of  this  roughness.  If  this  is  done,  the  answer 
to  the  problem  will  be  determined  in  a statistical  way.  This 
problem  was  studied  by  Patureau.  The  purpose  of  this  thesis 
is  to  test  the  above  theories  by  conducting  experiments  on  a 
laboratory  scale  model  and  comparing  these  results  with  field 
data.  Since  the  inertial  theory  was  unable  to  correlate  all 
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the  data,  another  theory  was  derived  from  the  results  of  known 
data. 

The  experimental  apparatus  and  procedure  will  be  studied 
In  detail  in  Section  2.  Since  the  results  of  Paturcau  are  ex- 
pressed in  terms  of  parameters  derived  from  the  power  spectrum 
of  the  surface,  this  section  also  described  the  apparatus  and 
procedure  used  to  determine  the  power  spectrum  of  the  experi- 
mental surface.  In  Section  3,  the  details  of  Patureau's  final 
size  theory  and  the  spreading  theory  derived  from  the  experi- 
mental results  will  be  discussed,  along  with  the  assumptions 
and  limitations.  Section  4 will  present  ihe  results  in  the 
form  of  correlations  between  theory  and  experiment,  an  explana- 
tion of  the  derivation  of  the  spreading  theory,  sources  of 
error,  and  the  conclusions  of  the  study. 


2.  EXPERIMENTS 


v In  planning  the  experiments,  the  characteristics  of  Arctic 
ice  were  first  studied.  Since  earlier  work  indicated  that  the 
spreading  was  gravity-inertia  dominated  it  was  decided  to  con- 
struct an  apparatus  that  could  dump  a known  volume  of  'oil'  at 
a fast  rate. 


To  obtain  dynamic  and  geometric  similarity,  certain  param- 
eters were  then  established  and  values  determined.  These  are: 


Re 

6 


(Ag<h~>)1/2<h~> 

V 


> 500 


<h  > 


1/80  above  ice 
1/20  below  ice 


v. 

N ■ number  of  pockets  >>  1 
A ■ density  ratio  <.2  (under  only) 
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With  the  exception  of  Re,  the  conditions  are  all  well  met.  In 
the  case  of  the  Reynolds  number,  a desire  to  use  fluids  that 
could  be  used  easily,  would  not  destroy  plexiglas  and  would 
provide  agreement  with  the  other  parameters  led  to  a relaxation 
of  the  restriction.  Water  for  the  over  ice  experiments  provides 
an  Re  of  a little  more  than  400  and  kerosene  for  the  under  ice 
case  provides  an  Re  of  about  100.  This  is  believed  to  provide 
sufficient  turbulence  to  satisfy  the  theory. 

A schematic  of  the  spreading  experiment  apparatus  is  shown 
in  Fig.  1.  This  figure  should  be  referred  to  if  any  questions 
arise . 

The  rough  surface  itself  was  built  by  epoxylng  various 
Sizes  of  rocks  to  a plywood  board  in  a random  fashion.  The 
surface  was  then  painted  white  with  black  concentric  circles. 
When  used  in  conjunction  with  a dyed  fluid,  the  spread  is  easily 
observed  and  photographed.  The  size  of  the  surface,  and  there- 
fore the  size  of  the  tank,  is  4*  x 4*.  It  was  felt  that  this 
size  tank  would  offer  a large  enough  surface  to  be  statistically 
significant  and  yet  small  enough  to  be  easily  handled. 

In  order  to  scale  the  laboratory  model  to  a real  spill,  one 
should  consider  that  a real  spill  would  be  on  the  order  of  10* 
cubic  feet  and  the  Arctic  roughness  height  is  on  the  order 
of  1'.  The  roughness  height  of  the  experimental  surface  is 
about  1/10  inch.  This  means  that  the  4'  x 4*  surface  cor- 
responds to  about  500  x 500  in  the  Arctic. 

The  tank  is  made  of  plexiglas  to  allow  the  spread  to  be 
easily  observed  and  photographed.  The  frame  was  designed  to 
allow  for  rotation  of  the  tank  when  conducting  under  ice  experi- 
ments. This  was  necessary  because  if  the  tank  were  filled  in 
a horizontal  position  the  pockets  in  the  surface  would  trap 
air  and  make  the  experimental  results  meaningless. 
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FIG.  1 SCHEMATIC  OF  SPREADING  APPARATUS  (OVER 


ICE  SPREAD) 


To  record  Che  progress  of  each  experiment  a Bolex  16mm 
camera  with  a 10mm  lens  and  a variable  time  base  was  used. 

Since  this  lens  needs  a focal  length  of  around  nine  feet  to 
get  a large  enough  field  to  cover  the  apparatus,  it  was  decided 
to  use  a mirror.  The  camera  was  then  run  at  a known  speed, 
usually  20  fps,  to  record  the  spread.  The  developed  film  was 
then  examined  frame  by  frame  to  determine  radius  versus  time 
and  the  final  size. 

To  measure  the  surface  a frame  was  made  from  1"  x 1"  angle 
iron  (see  Fig.  2)  and  marked  in  1/4"  increments.  The  elevation 
of  the  surface  was  then  measured  every  1/4"  with  an  accuracy  of 
1/32".  The  distance  below  a fixed  line,  the  frame,  was  re- 
corded by  hand  and  then  transferred  to  computer  cards.  The 
variance  and  the  power  spectrum  were  then  calculated  by  the 
program  shown  in  Fig.  3.  The  quarter  inch  spacings  were  chosen 
as  a compromise  between  accuracy  and  ease  of  measurement. 

To  test  the  isotropy  assumption,  traces  were  taken  in  per- 
pendicular directions  and  then  their  power  spectra  compared. 
This  is  explained  in  Section  4. 

The  frame  is  believed  to  affect  the  spectrum  only  in  the 
wavelength  range  greater  than  ten  inches.  If  this  is  true,  and 
there  are  no  indications  otherwise,  then  the  imperfections  in 
the  frame  will  have  no  effect  on  the  results. 

For  the  over  ice  experiments,  the  fluid  used  was  water 
dyed  with  food  coloring.  Soap  was  added  to  reduce  the  surface 
tension  and  the  'oil'  was  placed  in  the  cannister.  The  appara- 
tus was  then  leveled  very  carefully,  to  prevent  the  fluid  from 
running' off  to  one  side.  Problems  were  encountered  in  this 
irea.  Once  the  experiment  was  ready  to  start,  the  Camera  was 
started  and  the  stopper  pulled.  This  released  the  oil  at  a 
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FIG.  3 COMPUTER  PROGRAM  TC  DETERMINE  POWER 


SPECTRUM 


fairly  fast  rate.  In  order  to  determine  the  amount  released, 
the  fluid  level  in  the  cannister  was  measured  before  and  after 
the  release. 

The  underwater  experiments  were  considerably  more  difficult. 
The  experiment  had  to  be  set  up  in  such  a way  that  no  air  pockets 
remained  in  the  tank  and  no  kerosene  escaped  the  feed  lines  con- 
necting the  tank  to  the  cannister.  To  do  this,  water  was  first 
placed  in  the  cannister  to  a level  that  would  prevent  the  kero- 
sene from  getting  into  the  feed  lines.  Then,  enough  kerosene, 
which  was  dyed  using  enamel  paint,  is  added  to  fill  the  can- 
nister. The  rough  surface  was  then  placed  on  the  tank  and  the 
tank  was  almost  completely  filled  and  then  tilted  to  allow  the 
air  to  escape  from  the  pockets.  The  filling  of  the  tank  was 
completed  in  this  position.  When  all  the  air  escaped  out  the 
top  and  the  tank  was  full,  it  was  rotated  back  to  a horizontal 
position  and  leveled.  Since  there  was  sufficient  pressure  in 
the  tank  to  cause  the  surface  to  deflect  at  the  center,  water 
was  bled  off  in  an  attempt  to  ensure  that  the  surface  was  flat. 
Some  degree  of  error  was  probably  caused  by  this  deflection. 

Once  the  tank  was  leveled  and  the  pressures  equalized,  the  ex- 
periments proceeded  in  the  same  way  as  the  over  ice  experiments. 

To  check  for  the  presence  of  viscous  effects,  experiments 
were  run  with  glycerin-water  and  kerosene-oil  mixtures.  These 
experiments  are  described  in  the  Sub-Appendix. 

THEORY 

The  theory  that  governs  the  spreading  was  developed  as 
a result  of  data  from  the  experiments  described  in  the  last 
section  and  the  available  field  data.  The  steps  taken  in  de- 
riving the  theory  are  described  in  Section  4.  In  this  section 
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Che  theories  will  be  explained  in  simple  terms  and  the  analyti- 
cal expressions  derived. 

In  determining  what  will  happen  to  oil  spilled  in  the  Arctic, 
two  cases  must  be  considered.  In  the  case  of  a pipeline  break- 
age, the  oil  will  spread  over  ice  and  tundra.  In  the  case  of  a 
supertanker  accident,  the  oil  will  probably  spread  under  the  ice. 

This  happens  because  Arctic  ice  has  a specific  gravity  of  about 
0.8  and  North  Slope  crude  has  a specific  gravity  of  about  0.9. 
Therefore,  at  an  ice-water  interface  the  oil  will,  on  hydro- 
static considerations,  tend  to  spread  under  the  ice.  The  spread- 
ing in  both  cases  follows  the  same  theory  and  the  results  differ 
only  by  a factor  determined  by  the  density  ratio  of  water  and 
oil.  For  either  a pipeline  break  or  a supertanker  accident, 
the  maximum  rate  of  oil  released  will  be  about  1.5  x 10^  cubic 
feet  per  day. 

Patureau's  theory  states  that  pockets  in  the  ice  are  filled 
to  the  roughness  height,  <h“>,  and  it  was  shown  experimentally 
that  the  value  of  the  viscosity  is  not  important  (see  Sub-Appendix  B-l) . 
This  means  that  as  pockets  in  the  ice  are  filled  additional 
fluid  flows  over  the  filled  pockets  frictionlessly  and  the  drag 
can  act  only  at  the  edge  over  an  area  proportional  to  the  radius 
multiplied  by  the  roughness  height.  When  attempting  to  correlate 
the  laboratory  data  for  over  and  under  ice  and  the  field  data, 
it  was  found  that  the  data  would  not  scale  with  the  volume  re- 
leased. By  forming  dimensional  plots  described  in  the  next 
section,  it  was  shown  that  the  rate  of  spread  for  a given  sur- 
face could  be  scaled  using  the  rate  of  oil  release,  Q.  Since 
the  size  of  the  final  pool  depends  upon  various  statistical 
parameters,  one  can  see  that  for  a given  volume  and  surface 
there  vixl  be  an  average  pool  size  and  an  associated  variance. 

In  order  to  determine  the  accuracy  of  the  averages,  Patureau 
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developed  an  expression  relating  the  standard  deviation  of  the 
size  with  the  average  size. 


Figure  4 shows  a sketch  of  oil  being  released  on  a rough 
surface  with  a roughness  height  <h~>.  The  driving  force  is 
gravity,  which  causes  a pressure  at  the  outside  edge  of  the 
spill.  Thus, 


F ~ pA  ~ (pAgli)hr 
© 


Since 


h2.M>: 

rA 


the  gravitational  force  from  Eq.  (1)  can  be  expressed  as 


F * pAg 
8 r 


The  pressure  drop  caused  at  *he  outer  edge  as  the  oil  flows 


will  be 


dp  - Cf  pU 


■■■  ,:u.- 


where  Is  a friction  coefficient  of  order  one.  Since  the  area 
over  which  the  drag  acts  is  the  roughness  height,  <h~>,  times  r 
the  retarding  force  becomes 


Fr  ~ pU  r <h  > 


Replacing  U in  Eq.  (5)  with  r/t  and  setting  Eq.  (5)  equal  to 
Eq . (3)  gives 


pAg  -£212  - p ll  <h"> 
r t 


After  some  algebra,  this  becomes 


r - ( -S2_  ] t2/3 


This  can  be  non-dimensional ized  by  defining  a scaling  length 


r - <h  > 


and  a scaling  time 
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To  obtain  an  expression  for  the  final  size  and  the  associ- 
ated standard  deviation  additional  assumptions  must  be  made. 

The  first  of  these  assumptions  equates  the  average  of  a sample 
from  the  surface  with  the  ensemble  average  over  the  whole  sur- 
face. This  assumption  has  no  real  relevance  to  experimental 
procedure.  Next  it  is  assumed  that  the  elevation  of  points  on 
the  surface  are  normally  distributed  about  some  mean.  This 
assumption  is  in  fair  agreement  with  the  experimental  surface 
(and  good  agreement  with  the  real  ice  surface)  which  actually 
is  closer  to  a lag-normal  distribution.  Under  this  assumption, 
the  roughness  height,  <h  >,  and  the  variance  of  the  height,  c^, 
can  be  related  by  the  simple  equation: 


Thirdly,  it  is  assumed  that  pockets  in  the  ice  surface  are  filled 
to  the  roughness  height  (see  Fig.  5).  The  level  assumption  has 
no  real  justification  save  that  it  eases  the  calculations.  The 
accuracy  of  this  assumption  can  only  be  determined  by  experiment. 
The  fourth  main  assumption  is  that  the  number  of  pockets  filled 
is  always  large  enough  to  give  a statistically  meaningful  result. 
This  assumption  is  well  satisfied.  It  was  further  assumed  that 
the  pockets  assumed  a rectangular  shape  (square  in  the  case  of 
an  isotropic  surface).  This  was  done  strictly  for  computational 
reasons.  For  the  purpose  of  this  thesis  the  surface  will  be  as- 
sumed iso';ropic  and  the  spill  will  be  circular.  The  isotropy 
assumption  will  be  shown  to  be  a good  assumption  by  comparing 
the  power  spectra  along  perpendicular  paths  of  the  experimental 
surface.  Lastly,  it  was  assumed  that  a two-dimensional  power 
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spectrum  could  be  formed  as  the  product  of  two-dimensional 
spectra.  This  assumption  has  no  mathematical  justification, 
but  even  if  it  were  poor,  no  real  error  in  the  results  would 
appear.  The  procedure  for  determining  the  power  spectrum, 
which  is  done  by  computer,  is  explained  in  the  experimental 
section. 

The  final  results  of  the  theory  are  expressed  in  terms  of 
parameters  which  are  derived  from  the  statistics  of  the  ice 
surface.  The  first  of  these  is  the  variance  of  the  roughness, 
c^.  Next  are  two  parameters  k and  T derived  from  the  power 
spectrum  of  the  surface.  For  a one-dimensional  spectrum,  T 
is  the  wavelength  of  the  peak  of  the  spectrum  and  k is  deter- 
mined from  the  following  equation: 
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h , -2tt  . , , 2tt 

exp  ( - jj-j)  cosh  ( 


Tk 


Tk 


T72J 


(10) 


To  meet  the  requirements  of  the  theory,  the  product  kT  should 

2 2 

be  greater  than  50,  or,  in  the  two-dimensional  case  (kT  ) should 

^ 2 

be  greater  than  2500.  From  (kT“)  , an  additional  parameter,  0, 
can  be  determined  from  Fig.  6.  (See  eq  22) 

Following  the  assumpcion  that  the  oil  fills  the  pockets 
to  the  roughness  height  on  the  average,  we  get  the  equation 


<V>  - <S  > <h  > 

max 


or,  for  a normal  distribution 


<V> 

from  Eq.  (8) . 
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To  non-dimensionalize  the  results,  Patureau  Introduced  the  con- 
cept of  the  'most  probable  pocket'  size,^,  that  is,  the  most 
common  size  for  a pocket  on  the  ice.  Defining  the  surface  area 
of  a most  probable  pocket  as 


S - T2 
mp 


(13) 


and  the  volume  of  a most  probable  pocket  as 


mp 


T2  2<2V 


1/2 


(14) 


we  find  the  following  dimensionless  equation  from  Eq.  (12): 


<S  > 
max 


mp 


(±1 \ JL 

[T71) 


(15) 


To  provide  a measure  of  the  accuracy  of  the  above  equation, 
an  expression  for  the  variance  of  the  surface  area,  Os,  was 
derived.  This  expression  is  obtained  by  considering  the  pockets 
in  the  ice  to  be  filled  sequentially  and  independently.  For 
ease  of  understanding,  the  derivation  will  be  for  the  one- 
dimensional  case.  This  can  then  be  easily  extended  to  the  two- 
dimensional  case. 


Since  the  pockets  each  have  a volume  variance,  V , and 
a raP 

a length  variance  X and  since  they  are  filled  independently, 
the  total  variance  °V  will  be  expressed  as  the  sum  of  the 
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variances  of  Che  pockets  filled: 


Oy  - NOy  (16) 

mp 

The  total  length  of  the  average  one-dimensional  spill,  <X  >, 

luaX 

is  equal  to  the  length  of  a mean  pocket,  X^,  multiplied  by  the 
number  of  pockets  filled,  N. 


<X 

max 


> 


(17) 


°V 


(18) 


Using  Eq.  (16) 


a 


x 


(19) 


or 


In  dimensionless  form, 


this  becomes,  using  Eq.  (16): 


(19a) 
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In  the  two-dimensional  case,  one  finds 


o y2  /<S  >\ 

s [ max  ) 

s V 1 

mp  \ mp  / 


Where  the  proportionality  constant  can  be  defined  as 


" ( -J2E-  ) 
vmp  / 


with  0 as  shown  in  Fig.  6.  It  is  seen  that  0 depends  only  upon 


the  product  kT  . This  yields 


<S  > 
max 


which  is  the  result  presented  by  Patureau. 


4.  RESULTS 


Figures  7 and  8 show  the  power  spectra  for  surfaces  1 and 
2 respectively.  The  isotropy  assumption  can  be  checked  by  com- 
paring the  spectral  points  determined  from  paths  taken  per- 
pendicularly. The  two  directions  are  distinguished  by  the  use 
of  open  and  closed  symbols.  By  noting  that  there  is  little 
difference  between  the  two  directions,  one  can  see  that  the 
Isotropy  assumption  is  a good  one.  Also  shown  on  these  figures 


<7^,  k,  T and  0. 
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Figures  9 and  10  show  that  the  data  for  the  final  pool  size 

for  surfaces  1 and  2 respectively  is  in  good  agreement  with 

Patureau's  theory.  The  values  of  the  scaling  area  and  scaling 

volume  are  shown  on  the  figures,  where  S and  V are  as  calcu- 

tnp  mp 

lated  in  Eqs.  (13)  and  (14).  The  solid  lines  on  the  graphs  are 
the  average  of  the  final  size  from  Eq.  (15).  The  dashed  lines 
on  either  side  represent  one  standard  deviation  on  each  side  of 
the  average  and  were  determined  from  Eq.  (23).  On  these  graphs, 
the  open  symbols  represent  over  ice  runs  and  the  filled  symbols 
indicate  under  ice  runs.  It  appears  that  the  under  ice  pools 
are  slightly  smaller.  This  could  be  within  experimental  error 
or  could  be  caused  by  a slight  deflection  of  the  rough  surface 
caused  by  water  pressure.  Water  was  bled  from  the  tank,  in  an 
attempt  to  remove  the  deflections,  but  some  error  is  anticipated. 

Using  the  inertial  theory  to  correlate  the  laboratory  data 
led  to  two  distinct  groups,  over  and  under  ice,  with  intercepts 
differing  by  about  a factor  of  two.  Additionally,  the  over  ice 
data  did  not  correlate  well  with  the  data  from  Glaeser  and  Vance. 
In  order  to  determine  the  proper  scaling  laws,  a series  of  dimen- 
sional graphs  were  made.  Since  the  radius  in  the  laboratory  data 

appeared  to  grow  approximately  as  the  time  to  the  one-half  power, 

1/2 

it  was  decided  to  calculate  r/t  for  all  laboratory  and  field 

data.  These  values  were  then  plotted  against  volume  released, 

release  rate,  gravity,  and  roughness  height.  Figure  12  shows 
1/2 

r/t  as  a function  of  volume  on  the  left  side  and  release 

1/2 

rate  on  the  right.  It  can  be  seen  that  r/t  cannot  be  de- 
l/2 

scribed  as  a function  of  volume  but  that  r/t  does  vary  ap- 
proximately as  the  release  rate  to  the  one-third  power.  Once 
a rough  power  law  for  release  rate  was  known,  additional  correla- 
tions were  tried.  Using  the  results  of  these  crossplots  and  the 
viscosity  experiments,  the  theory  described  in  Section  III  was 
derived. 


FIG.  10  DIMENSIONLESS  VOLUME  RELEASED  vs  SURFACE 


AREA  FOR  SURFACE  2 
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Figure  13  shows  the  dimensionless  radius  versus  dimension- 
less time  for  both  over  and  under  ice  cases  using  the  results 
of  Section  3.  Also  shown  on  this  graph  is  field  data  from 
Glaeser  and  Vance ^ and  McMinn^.  The  raw  laboratory  data 
that  make  up  this  figure  is  given  along  with  the  dimensionless 
values  in  Sub-Appendix  B-3.  The  equation  that  best  describes  the 
line  in  Fig.  13  is: 

r - . 25 ( AgQ2 / <h”> ) 1 7 6 t2/3 

The  under  ice  data  shown  in  Fig.  13  is  plotted  after  allow- 
ing for  a starting  time.  It  was  found  that  in  the  under  ice 
case  there  was  a significant  lag  between  the  time  the  stopper 
was  pulled  and  the  actual  beginning  of  the  spread.  This  was 
due  to  the  time  for  the  kerosene  to  flu»t-  up  through  the  water 
to  the  surface.  This  starting  time  was  found  to  be  about  1.5 
seconds  (see  Sub-Appendix  B-2). 

For  the  field  data,  the  roughness  heights  are  not  known. 
Thus,  the  field  data  was  fit  to  the  solid  line  in  Fig.  13  by 
adjusting  the  value  of  the  roughness  height.  In  some  cases, 
this  led  to  a value  of  the  roughness  height  that  appears  either 
too  large  or  too  small. 

However,  the  experiments  were  run  in  small  areas  with  dif- 
ferent kinds  of  surfaces  (snow  covered,  windswept  ice,  etc.), 
so  effective  roughness  heights  for  these  experiments  could  vary 
greatly. 

To  expand  the  results  of  the  laboratory  experiments,  one 
must  know  the  characteristics,  particularly  the  roughness  height, 
of  the  Arctic  ice.  From  laser  profilometer  traces  taken  in  areas 
of  high  spill  potential  it  has  been  determined  that  the  rough- 
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FIG.  13  DIMENSIONLESS  RADIUS  vs  TIME  FOR  LABORATORY  AND  FIELD  DATA 
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ness  height  varies  greatly  but  is  seldom  less  than  three  centi- 


meters or  greater  than  60  centimeters.  The  size  of  a spill  would 

.4 


be  on  the  order  of  10  cubic  meters.  Figure  14  shows  how  the 
final  radius  of  a spill  varies  with  the  roughness  height  and 
the  total  volume  released.  For  an  average  volume  and  roughness 
height  the  final  radius  will  be  on  the  order  of  300  meters. 


110 


SUB-APPENDIX  B-l 


Since  the  Reynolds  Numbers  of  the  under  ice  and  over  ice 
experiments  differ  by  a factor  of  about  five  and  the  starting  time 
of  the  under  ice  and  over  ice  data  differ  by  about  60%,  it  was 
decided  to  test  for  viscosity  effects.  This  was  done  by  Increasing 
the  viscosity  of  the  over  ice  fluid  (soap  and  water)  by  the  addition 
of  glycerin  and  by  increasing  that  of  the  under  ice  fluid  (kerosene) 
with  SAE  30  wt.  oil.  The  viscosity  of  the  new  fluids  was  measured  on 
a Tag  Saybold  viscometer.  The  figures  for  the  viscosities  of  the 
various  experiments  are  shown  in  the  table  below.  Changing  the 
viscosity  of  the  water  by  a factor  of  five  and  the  viscosity  of 
kerosene  by  factors  of  two  and  four  had  no  real  effect  on  the 
fluid  spread. 

The  Reynolds  number  below  is  defined  as: 


npr  z 

- r O 0 

Re  . 

For  the  over  ice  experiments  with  just  water,  the  Reynolds  number 
was  about  400  and  for  the  kerosene  only  under  ice  experiments,  the 
Reynolds  number  was  about  80. 
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SUB-APPENDIX  B-2 

In  Che  case  where  Che  kerosene  was  released  under  Che  rough 
surface,  Chere  was  a delay  becveen  Che  Cime  Che  scopper  was  pulled 
and  Che  Cime  ChaC  Che  fluid  scarced  to  spread.  This  'starting 
time1  was  found  by  extrapolating  Che  radius  versus  Cime  curve  Co 
zero  radius.  Since  Che  radius  varies  approximately  as  Che  square 
root  of  the  elapsed  cime,  a graph  was  constructed  using  Che  radius 
as  the  ordinate  and  Che  cime  squared  as  the  Absiosa  (see  Fig.  15). 
By  extrapolating,  the  starting  Cime  was  found  Co  be  about  1.5 
seconds . 
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SUB-APPENDIX  B-3 


Below  la  the  raw  data  for  the  laboratory  experiments  along  with 
the  non-dimensionalized  values:  (in  cgs  units) 


Run  1 

Over 

Surf.  1^ 

Run  2^ 

Over 

Surf.  1 

Q - 88 

V 

- 90 

r - 16*5 

max 

Q - 270 

V - 218 

r - 24 

max  — 

r 

f 

t 

t 

r ? 

t 

E 

10.5 

121 

.70 

2550 

12  138 

.5 

2595 

12 

138 

1.1 

4005 

13  150 

.65 

3375 

14 

161 

1.5 

5460 

15  173 

1.0 

5190 

17  196 

1.45 

7525 

19  219 

2.0 

10.380 

Run  3_ 

Over 

Surf.  1^ 

Run  4^ 

Over. 

Surf.  1. 

Q - 190 

V 

- 99 

r - 18 

max  — 

Q - 123 

V - 104 

r - 17 

max  — 

r 

r 

t 

t 

r ? 

t 

t 

11 

127 

.65 

2840 

10  116 

.60 

2110 

13 

150 

.80 

3490 

13  150 

.85 

2980 

15 

173 

1.3 

5680 

14  162 

1.05 

3690 

Run  5 Under 

Q - 22.5  V - 144 
r r t 

Surf.  1^ 

r - 18 

max  — 

t 

Run  6 

Q - 21^ 
•* 

x r 

Under 
V - 126 
t 

Surf.  jL 

r - 18.5 

max  

** 

t 

10.5  121 

3.1 

2890 

9.3 

110 

3.0 

2800 

12  138 

4.3 

4010 

12 

138 

4.3 

4015 

13  150 

5.5 

5130 

13 

150 

5.1 

4765 

15  172 

7.0 

6525 

15 

172 

8.5 

7940 

Preceding  page  blank 
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SUB-APPENDIX  B-3 
(continued) 


Run 

Q - 

i 7 Under 

21.9  V - 153 

Surf.  2 
r - 19.5 

Run  8 Under 

Q - 20.6  V - 99 

Surf.  2^ 
r • 16 

r 

? t 

BAX 

t 

r 

? t 

BBX  — 

t 

10 

105  3.7 

3030 

9.0 

95  2.8 

2210 

12 

126  5.6 

4585 

12 

126  5.3 

4185 

15 

158  8.8 

7207 

15 

158  8.8 

6952 

Run 

9 Under 

Surf.  2^ 

Run  14  Over 

Surf.  2 

Q - 

22.5  V - 117 

r *16 

max  — 

A* 

Q - 

180  V - 126 

r - 20 

max  — 

A* 

r 

? t 

t 

r 

r t 

t 

9 

95  3.8 

3155 

12 

126  .6 

2230 

12 

126  4.2 

3490 

14 

147  1.0 

3710 

14 

147  5.9 

4900 

15 

158  1.25 

4640 

Run 

15  Over 

Surf.  2_ 

Run 

. 14  Over 

Surf.  2_ 

Q - , 

180  V - 135 

r - 19 

BAX  

Q - 

180  V - 126 

r - 19 

max  — 

r 

? t 

l 

r 

i t 

E 

10 

105  .4 

1480 

10 

105  .4 

1460 

13 

137  .75 

3180 

12 

126  .55 

2040 

15 

158  1.1 

4660 

14 

147  1.05 

3900 

Run 

17  Under 

Surf.  2 

Run 

18  Under 

Surf.  2 

Q - , 

14.4  V - 126 

a# 

r - 

BAX 

A* 

Q - 

20.0  V - 180 

A* 

r ■ 

max 

A# 

r 

r t 

t 

r 

r t 

t 

10 

105  3.5 

3600 

10 

105  2.8 

3540 

12 

126  4.9 

5080 

12 

126  3.8 

4810 

15 

158  6.5 

6730 

15 

158  5.8 

7340 

i 

3 
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APPENDIX  C 
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STATISTICAL  APPROACH 
EXTENT  OF  AN  OIL  SPILL 


FOR  DETERMINING  THE 
OVER  A ROUGH  SURFACE 


1.  INTRODUCTION 

As  one  seriously  thinks  of  using  the  crude  oil  discovered 
in  the  Arctic  Slope,  no  natter  which  means  of  transportation 
will  be  adopted,  it  appears  that  some  understanding  of  how  oil 
would  spread  on  the  ice  is  of  major  importance. 

This  problem  cannot  be  dealt  with  as  a classical  oil  spill 
on  a plane  surface  because  the  characteristic  feature  of  an 
Arctic  Ocean  ice  surface  is  its  roughness.  Moreover,  field  ob- 
servations show  that  such  a surface  can  be  regarded  as  formed 
of  'pockets"  of  different  sizes,  some  of  them  being  intercon- 
nected. One  can  easily  see  that  the  oil  will  stop  spreading 
abruptly  after  having  filled  a certain  number  of  pockets.  The 
goal  of  this  work  is  to  determine  the  final  size  of  the  pool 
as  a function  of  the  amount  of  oil  released. 

It  is  clear  that  the  geometry  of  the  ' pockets"  previously 
introduced  is  entirely  related  to  the  local  variations  in  the 
ice  surface  altitude;  it  can  therefore  be  determined  by  func- 
tions h(x,y)  defined  in  an  arbitrary  plane  coordinate  system 
xy  as  the  elevation  of  the  ice  surface  above  a reference  level; 
most  often  in  the  rest  of  this  work,  these  functions  will  be 
referred  to  as  roughness  height.  A simple  way  to  record  informa- 
tion about  these  variations  in  the  elevation  is  to  take  airborne 
laser  profilometer  traces^.  Such  data  were  obtained  by  the  U.S. 

A profiler  operating  from  a plane  records  the  phase  delay  be- 
tween transmitted  to  the  ground  and  reflected  laser  beam  light, 
thus  measuring  the  altitude  variations  along  straight  lines  and 
yielding  a random  output  signal. 
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Navy  in  the  Beaufort  and  Chuckchi  Seas.  There  is  no  doubt, 
when  examining  any  of  these  traces,  that  the  major  character 
of  the  ice  roughness  is  its  randomness  (Ref.  [3]).  Releasing 
the  same  amount  of  oil  V at  different  locations  on  a randomly 
rough  ice  surface  would  result  in  different  values  for  the 
measurement  of  the  oil  pool  final  size.  Clearly,  we  must  then 
formulate  our  problem  in  a statistical  way,  namely:  the  extent 

to  which  oil  would  spread  has  to  be  determined  in  terms  of  a 
certain  number  of  statistical  parameters  characterizing  the 
ice  surface  where  the  spill  occurs.  An  analytical  treatment 
of  the  profilometer  traces  can  provide  us  with  h^,  the  mean 
value  of_the  roughness  height  h,  a.,  the  variance  of  h,  equal 
to  (h  - h)  , and  the  auto-correlation  function  of  h.  We  will 
see  in  the  course  of  this  work  that  these  quantities  are  suf  - 
ficient to  give  an  answer  to  our  problem,  provided  certain 
simplifying  assumptions  are  made. 

With  no  more  field  data  than  those  just  introduced,  it  is 
unlikely,  however,  that  we  should  end  up  with  an  elaborate  and 
complete  description  of  the  oil  spreading  phenomenon.  This 
paper  will  only  aim  at  determining  the  first  two  statistical 
moments  of  the  oil  pool  final  size  - i.e.,  its  average  value 
and  its  variance  - both  in  cases  when  the  spreading  follows  a 
privileged  direction  (so  called  "one-dimensional")  and  when  it 
doesn't  (two-dimensional). 

Section  2 will  introduce  and  discuss  the  basic  assumptions 
necessary  to  set  up  a simplified  model  of  the  phenomenon,  where- 
as Section  3 will  give  a general  mathematical  expression  for 
these  quantities  in  terms  of  the  field  data  available,  namely: 
h,  0^  and  the  auto-correlation  of  h. 


The  bar  denotes  an  average  over  a sample. 
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In  Section  4,  we  will  be  concerned  with  applying  the  pre- 
vious results  to  certain  ice  surfaces  where  field  data  Is  avail- 
able, thereby  involving  the  modelling  of  the  auto-correlation 
function  of  h or,  more  exactly  of  its  cosine  Fourier  transform. 
Finally,  in  those  specific  cases  examined  in  Section  IV,  a 
simple  power  law  will  be  derived,  relating  the  first  two  stat- 
istical moments  of  the  final  size  of  the  oil  pool  In  « dimension- 
less form. 

2.  BASIC  ASSUMPTIONS 

If  the  ice  surface  where  the  spill  occurs  is  being  con- 
sidered as  randomly  rough,  it  is  clear  that  we  cannot  have  a 
complete  definition  of  the  function  h relevant  to  that  specific 
spill,  which  is  what  would  be  necessary  to  yield  a single,  defi- 
nite value  for  the  oil  pool  size  in  that  case. 

A way  of  overcoming  this  difficulty  is  to  retain  a certain 
number  of  statistical  properties  of  h met  by  different  possible 
"realizations"^  of  the  ice  surface.  These  properties  are  then 
regarded  as  "ensemble  averaged"  over  all  those  realizations  and 
generally  consist  of  a certain  number  of  locally  defined  stat- 
istical moments  and  correlation  functions  of  h.  These  quantities 
provide  a certain  amount  of  information  common  to  any  ice  surface 
realization  that  belongs  to  the  ensemble.  Obviously,  the  smaller 
this  number,  the  more  random  is  the  surface  and,  when  the  num- 
ber of  such  properties  goes  to  infinity,  the  ice  surface  tends 
to  be  perfectly  defined  and  therefore  no  longer  random.  In 
other  words,  since  we  are  talking  about  ensemble  averages,  all 


A "realization"  of  the  ice  surface  verifying  certain  statistical 
conditions  is  an  actual,  observed  surface  such  that  it  meets  the 
previous  conditions. 
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(he  realizations  are  becoming  identical.  The  idea  is  thus  to 
haye  locally  defined  functions  that  give  statistical  informa- 
tion on  a whole  ensemble  of  possible  realizations  of  the  ice 
surface. 

In  the  present  case,  the  only  statistical  data  we  have  is 
_ _ 2 

h,  < ■ (h  - h)  and  h(P)h(P')  where  the  bar  denotes  an  average 
pver  a specific  "realization"  of  the  ice  surface  and  where  P 
and  P*  denote  any  two  different  points  of  that  surface.  Since 
these  data  are  "space  averaged"  and  therefore  not  locally  de- 
fined, we  are  lead  to  make  an  ergodic  hypothesis,  assimilating 
77  and  <.  ,>^  averages.  Our  random  surface  is  now  characterized 
by  <h>,  0^  ■ <[h  - h ]^>  and  h(P)h(P')  but  clearly,  this  is 
not  sufficient  since  nothing  is  said  about  the  aspect  of  the 
surface. 

There  appears  our  second  assumption,  related  to  the  stat- 
istics of  the  roughness  height  h:  we  will  suppose  h is  normally 
distributed.  If  6(h)  denotes  this  distribution,  we  have  im- 
mediately: 


6(h) 


/2ir  o' 


exp 


(h  - <h>) 

2a. 


(1) 


6(h)dh  represents  the  fraction  of  time  that  h will  be  found  be- 
tween h and  h + dh.  It  does  not  seem  possible  to  prove  this 
assumption  by  any  physical  argument  since  the  natural  processes 
that  make  the  ice  grow  are  extremely  intricate  and  not  thoroughly 
known.  Therefore,  we  have  to  verify  in  each  case  that  this  is 
not  too-  far  from  the  observed  distribution;  this  was  done  as 


ensemble  and  space  averages  are  equal. 
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an  example  on  Fig.  3 where  the  Gaussian  curve  J 6(x)dT  is 

•60 

plotted  against  h and  is  compared  to  corresponding  data  points 
measured  on  a real  surface.  It  appears  that  the  Gaussian  ap- 
proximation is  in  relatively  good  agreement  with  observations. 

The  third  assumption  concerns  the  way  the  oil  would  fill 
the  "pockets"  of  which  the  ice  surface  is  formed  and  more  pre- 
cisely, up  to  what  level.  Since  the  ice  relief  is  random,  the 
answer  to  this  question  is  unlikely  to  be  unique.  Most  probably, 
a large  number  of  different  levels  would  appear,  corresponding 
to  more  or  less  "open'1  or  "closed"  pockets.  However,  to  keep 
the  problem  feasible,  we  have  to  account  only  for  an  average 
uniform  level  of  the  oil  pool,  which  is  another  way  to  say  that 
all  the  "pockets"  are  interconnected.  Itoreover , it  seems  a 
reasonable  assumption  - at  least  most  simplifying  - to  say  that 
this  uniform  level  is  the  same  as  that  of  <h>.  No  observations 
exist  so  far  to  confirm  such  a point  and  the  reason  why  we 
chose  this  level  is  that  any  other  one  would  have  been  even 
more  arbitrary  and  would  have  lead  to  more  intricate  mathemati- 
cal derivations. 

The  fourth  hypothesis  reads  as  follows : we  will  always 

assume  throughout  this  study  that  the  number  of  "pockets"  in 
volved  in  an  oil  slick  is  very  large;  in  other  words,  the 
ratio  of  the  volume  V released  to  the  mean  volume  of  a pocket 
is  large;  this  will  allow  us  to  give  a statistical  answer 
for  the  final  size  of  the  oil  pool.  This  assumption  may  be 
considered  as  valid  in  so  far  V is  generally  large  (10,000  to 
50,000  m^  for  the  crash  of  a supertanker)  and  also  the  ice 
roughness  is  such  that  a big  pocket  contains  many  small  ones 
and  so  on;  therefore,  even  if  the  slick  occurs  in  a big  one, 
the  numbers  of  pockets  involved  will  always  be  large. 
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These  four  assumptions  obviously  cannot  lead  to  a sophisti- 
cated model  since  they  mainly  reflect  our  lack  of  information 
about  ice  topography  in  the  Arctic  Slope.  Moreover,  the  general 
mathematical  results  obtained  from  this  model  won't  be  the 
unique  approximate  solution  to  the  oil  spreading  problem  since, 
for  instance,  another  frequency  distribution  for  h could  have 
been  selected  - log-normal  particularly  -,  or  another  height 
of  the  oil  pool  level  could  have  been  chosen.  Nevertheless, 
the  simple  model  we  can  thus  build  provides  us  with  a first 
understanding  of  how  oil  spreads  over  a statistical  surface. 


3.  MATHEMATICAL  ANALYSIS 

The  first  four  subsections  will  be  concerned  with  a two- 
dimensional  study  whereas  subsection  3.5  will  briefly  state 
the  results  for  a one-dimensional  case. 

3.1  Problem 

We  suppose  no  privileged  direction  in  the  ice  sur- 
face, at  least  for  the  scale  of  the  roughness  we  are  in- 
terested in^.  Therefore,  the  ice  topography  is  determined 
by  a two-dimensional  function  h(x,y)  describing  the  vari- 
ations of  the  surface  altitude  with  respect  to  a reference 
plane.  The  origin  of  the  coordinates  x and  y is  taken  at 
the  point  of  oil  release  and  the  reference  level  is  chosen 
so  that: 

<h(x,y)>  - 0 (2) 

It  has  been  shown  that  big  pressure  ridges  do  appear  with 
certain  definite  orientations. 
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There  again,  <h(x,y)>  denotes  the  ensemble  average  of 
h(x,y)  where  the  averaging  ensemble  consists  of  the  many 
lea  surface  realizations  that  can  be  observed  in  the  region 
one  considers.  From  our  hypothesis,  <h(x,y)>  does  not  de- 
pend on  x and  y and  therefore  has  to  be  identified  to  the 
space  average  h(x,y)  of  h(x,y)  over  the  realizations  of 
the  ice  surface  themselves. 

2 

By  the  same  token,  <?h  ■ <h  (x,y)>  is  the  variance  of 
h(x,y)  calculated  over  the  whole  ensemble,  lienee,  the 
frequency  distribution  of  h(x,y),  assumed  to  be  Gaussian, 
is  independent  of  any  direction  and  is  given  by  Eq.  (1). 

A proper  investigation  of  the  actual  ice  topography 
f 

can  provide  us  with  a two-dimensional  auto-correlation 
function  defined  as: 


'f(t.n)  - <h(x  + £,  y + n)  h (x,y)>  (3) 


We  can  thus  pose  the  two-dimensional  problem  in  the 

following  terms:  in  a region  of  the  Arctic  Ocean  where 

2 

we  know  0^  ■ <h  (x,y)>  and  4'(C,n)t  a volume  V of  oil  is 
released  at  point  x ■ y • 0.  What  will  be  the  total  area 
covered  with  oil  and  also,  what  will  be  the  final  geo- 
metrical aspect  of  the  pool? 

The  general  mathematical  answer  will  consist  of  the 
first  two  statistical  moments  of  the  oil  pool  size  in 


+ 

By  building  a two-dimensional  grid  of  laser  profilometer 
traces. 
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terms  of  V,  c.„  and  f(C>n)i  along  with  a dissymmetry  factor 
expressing  the  most  probable  ratio  of  the  longest  dimen- 
sion of  the  pool  to  the  smallest  one. 


3.2  Dissymmetry  of  the  Final  Oil  Pool 


As  we  do  not  allow  any  dissymmetry  in  the  frequency 
distribution  of  h(x,y),  the  final  shape  of  the  slick  will 
be  influenced  by  the  aspect  of  the  "pockets"  themselves 
in  the  plane  xy.  In  general,  the  ice  surface  is  not  iso- 
tropically rough.  It  follows  that  if  Xx  and  Xy  are  the 
ensemble  averages  of  a pocket  dimensions  in  the  x and  y 
directions  respectively,  the  ratio  y of  Xx  to  Xy  can  be 
called  "dissymmetry  factor"  since  it  provides  a rough 
idea  about  the  final  pool  aspect. 

From  Sub-Appendix  C-l,  it  appears  that  y is  given  by: 


(0,0) 


Now,  let  us  imagine  how  oil  would  spread  from  point 
of  release.  Since  the  volume  V is  supposed  to  be  very 
large,  the  shape  of  the  pockets  is  not  important  at  the 
beginning  and  the  slick  would  propagate  circularly.  But, 
after  this  initial  phase,  one  can  assume  that  the  geometri 
cal  aspect  of  the  pockets  would  influence  the  final  shape 
of  the  oil  slick. 

Considering  a laser  profllometer  trace  that  passes 
through  the  origin  of  the  coordinates  x and  y,  it  is  clear 
from  our  ergodic  hypothesis  that  the  statistics  of  the 
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roughness  height  remain  the  same  all  along  the  trace.  It 
follows  that  the  origin  of  the  coordinates  is  a geometric 
center  of  symmetry  for  the  pool  shape.  Moreover,  in  order 
to  simplify  further  calculations,  we  assume  that  the  re- 
sulting oil  pool  is  rectangular,  xmgx  and  ymax  being  the 
x-direction  and  y-direction  dimensions  respectively. 

Therefore,  it  follows  that: 


* y i unity,  in  general 


Obviously,  when  the  ice  surface  is  truly  isotropic,  y is 
unity  and  the  final  pool  is  a square  in  our  assumption. 


3 . 3 First  Two  Statistical  Moments  of  the 
Oil  Pool  Final  Size 

As  was  suggested  at  the  end  of  the  last  sub-section, 
the  oil  pool  is  likely  to  be  symnetric  with  respect  to  the 
point  of  release  i.e.»  the  origin  of  the  coordinates. 
Hence,  defining  h+(x,y)  and  h (x,y)  such  that  h - h+  + h 


h~(x,y)  - 0 


when  h(x,y)  £ 0 


h (x,y)  ■ h(x,y)  when  h(x,y)  < 0 


we  can  introduce  v(x,y)  as  the  volume  held  by  the  ice  sur- 
face - under  the  third  assumption  of  Section  3 - within  a 


rectangle  of  dimensions  x and  y.  Then: 


n/2 

h“(x,y)  dx  dy  (6) 

_ 


This  is  a mathematic  definition  and,  evidently,  it  does 
not  account  for  the  fact  that  the  oil-air  interface  must 
be  horizontal  when  edge  pockets  are  not  completely  filled 
up  to  the  average  uniform  level;  this  fact  is  illustrated 
on  the  sketch,  Fig.  1 for  the  one-dimensional  spreading. 

It  follows  that  the  mathematical  extent  of  the  oil  is  only 
slightly  different  from  the  real  one:  at  most  within  a 

pocket.  Since  we  made  the  assumption  of  a large  number 
of  pockets,  we  can  define  x and  y as  the  x-direction 
and  y-direction  maximum  extents  of  the  final  pool  when  a 
volume  V of  oil  is  being  released  over  a specific  realiza- 
tion of  the  ice  surface;  namely: 

h”(x,y)  dx  dy  (7) 

Or,  introducing  the  area  SMx  defined  as: 


x /2 
max 


W 


max  'raa: 

-J  12  -y  /2 

max  •'max 


3 


Relation  (7)  becomes: 


■ - I h“(x,y)  dx  dy 


(9) 


max 


As  the  function  h is  everywhere  continuous,  we  can 
also  write  relation  (9)  as: 


V - - I h~(x,y)  dx  dy  - [S^  - <smax>lh  (y,p)  (10) 


<S  > 
max 


where  point  (y,p)  belongs  to  the  surface  (R  - <S  >) ; 
v max  max 


h (y,p)  is  exactly  the  average  of  h taken  over  the  sur- 


f*“  (!w  - * W>  • 


We  now  assume  that: 


h“(y,p)  • <h“(5,n)> 


(ID 


This  cannot  be  justified  in  general.  Therefore,  Eq.  (11) 
is  more  valid  as  this  surface  becomes  large  since,  when 
the  averaging  sample  gets  large  enough,  the  average  of  h 


approaches  h~(5,n),  hence  <h  (C *n) > from  our  ergodic  hy- 
pothesis. An  "a  posteriori"  verification  will  thus  be 
possible. 


Then,  taking  the  ensemble  average  of  (10)  and  making 
use  of  the  previous  approximation,  we  obtain  the  first 
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statistical  moment  of  the  pool  size,  namely: 


<S  > 
max 


V 

<h"(£,n)> 


(12) 


To  determine  the  second  moment,  let  us  first  notice 

that  o , the  variance  of  S and  a the  variance  of 
s max  v 

v(<x  >,  <y  >)  are  related.  By  definition: 
max  max  1 


a 

e 


<(S, 


max 


<S  >)2> 
max 


(13) 


C ■ < |v(<x  >,<7  >)  - <\>(<X  >,<y  >)>|2>  (14) 

V r ' max  * *max  7 v max  * ymax  7 J v 7 


From  relations  (10)  and  (11),  it  comes  out  that: 


J 


h (x,y)  - V 


<S  > 
max 


<(S  - <S  »2>  <h~(£,n)2  (15) 

max  max 


Or: 


o <h“(Cln)2> 


(16) 


r 


= i 


t- 

i 

i 


i 

t 


Therefore,  computing  in  Sub-Appendix  C-2,  we  end  up  with: 

<x  >//2  <y  >//2 

max  max  r.  . -1 

a - 2~  [ [ - K 

<h"(S,nr>  J J L ft 


[¥  ■ n]  ■ s] 


dn 


(17) 


From  all  previous  calculations,  especially  relations 

(12)  and  (17),  it  appears  that  knowing  <h  (£,n)>  and 

¥ (£»Tl)  enables  us  to  derive  the  average  value  <S  > and 

1/2  max 

the  standard  deviation  0 of  the  total  area  covered 

8 

with  oil,  in  terms  of  the  amount  V released.  Furthermore, 
an  idea  about  the  dissymmetry  in  the  aspect  of  the  pool 
can  be  gotten  from  coefficient  y.  This  will  constitute 
our  answer  to  the  two-dimensional  problem  posed  in  Sub- 
section 3.1,  provided  we  can  relate  <h  (£,n)>  and  f (£,n) 
to  the  available  statistical  data  i.e.,  0^  and  ^(^.n). 

This  is  the  purpose  of  the  next  sub-section. 


3.4  Determination  of  h (C,n)  and  (5,n) 

He  now  go  into  more  detail  about  the  distribution 
function  5(h)  that  was  first  mentioned  in  the  Introduction. 
This  function  tells  us  that  the  probability  h will  lie 
between  y and  y + dy  is  5(y)dy  and  therefore: 


(18) 
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Assuming,  as  In  our  second  hypothesis,  that  this  fre- 
quency distribution  is  Gaussian  with  mean  value  0 implies 
that  h will  be  negative  as  many  times  as  it  is  positive. 
Since  y is  the  value  taken  by  a stochastic  function  h,  it 
is  also  a stochastic  variable.  Therefore,  as  we  know  its 
frequency  distribution  function  6(y),  it  is  clear  that  the 
average  of  a stochastic  quantity  Q uniquely  related  to  y 
is  expressed  by: 


<Q> 


/■ 


Q(Y><$(Y)d 


(19) 


If  now  we  identify  Q to  the  value  of  h (x,y),  it  fol- 
lows from  relations  (5) : 


Q ■ y when  y 4 Q 

Q ■ 0 when  y > 0 (20) 


Then,  substituting  into  Eq.  (19)  for  6(y)  and  Q(y) : 


Or,  performing  the  integration: 


<h"(C,n)> 


(21) 


(22) 
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Nov,  to  determine  we  have  to  extend  what  we 

previously  did  for  a one-dimensional  distribution  of  h to 
H joint  two-dimensional  frequency  distribution.  We  end  up 
vith  the  following  relation  (see  Sub-Appendix  C-3  for  details) : 


1/2 


2 1 

(u  - j)  erf  (u)  exp 


where  we  have  set  for  convenience: 


J 

f 


V - ¥(0,0) 

T - ¥<5,n)  (24) 


It  appears  (see  relation  (2)  in  Sub-Appendix  C-4)  that 

¥” (0,0)  - y ¥(0,0)  which  is  consistent  with  the  physical 

conditions  at  ¥ - ¥ namely: 

o 


¥(0,0)  - <h2>  - <(h+  + h-)2> 


+ 2<h+  h“>  - 2<h"  > 
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> + <h 


2 

> 


(25) 
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When  assuming  the  statistics  of  h and  h+  are  the  same  as 
suggested  by  a normally  distributed  h.  The  product  h+h 
must  always  be  zero  as  can  be  seen  by  Eq.  (5). 

Furthermore,  when  C or  n tends  to  infinity,  ^ (£,n) 

°h 

tends  to  , in  compliance  with  the  fact  that  h (x,y) 
and  h (x  + 5,y  + n)  are  then  no  longer  correlated.  There- 
fore, according  to  relation  (22): 


<h“(x,y)h~(x  + ,y  + )>  - <h“(x,y)><h~(x  + £,y  + n)> 

- <h~(5,n)>2  - ah/2TT  (26) 


We  have  therefore  provided  a means  of  deriving  the  statisti- 
cal quantities  we  need  in  Sub-section  3.3,  i.e.,  ¥ (£,n) 

and  <h  (g,n)>,  from  the  field  data  generally  available, 

2 

namely:  0^  ■ (h  - h)  , the  variance  of  the  roughness  height 

end  f(£,n),  its  auto-correlation  function. 

3.5  One-Dimensional  Case 

We  can  imagine  that  in  a particular  region,  the  statistics 
of  the  ice  roughness  are  almost  one-dimensional;  that  is  to 
say,  the  averaging  ensemble  consists  here  of  traces  parallel 
to  this  privileged  direction  we  denote  by  x.  Therefore,  the 
ensemble  <h>  and  auto-correlation  function  are  only  func- 
tions of  x.  And,  if  h(x)  denotes  the  elevation  (positive  or 
negative)  above  the  mean  value  of  h taken  along  the  trace 
(see  Fig.  1),  we  have  from  our  ergodic  hypothesis: 
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<h(x)>  - h(x)  " 0 

<h2(x)>  - h2(x)  - o variance  of  h 


and,  ¥(x)  being  the  auto-correlation  function  of  h(x): 


t|»(T)  - <h(x)h(x  + T)>  - h(x)h(x  + f) 


From  the  previous  considerations,  it  follows  that  any 
actual  surface  profile  parallel  to  x will  be  described  by 
the  statistical  properties  (27)  and  (28). 

Now,  suppose  an  oil  spill  occurs  in  such  a zone,  a 
volume  V being  released  per  unit  width  perpendicular  to  x; 
if  the  spill  is  wide  enough,  we  can  regard  the  spreading 
of  oil  as  one-dimensional  in  the  direction  of  x.  This  is 
not  completely  unrealistic  since  narrow  water  leads  bounded 
by  large  floes  are  fairly  frequent  in  the  Arctic  Slope. 

Under  these  conditions,  a large  oil  spill  would  first  spread 
rapidly  on  water,  thus  creating  this  wide  spill  (V  per  unit 
width)  that  would  further  spread  over  ice  almost  in  a single 
direction  perpendicular  to  the  lead. 

We  can  then  put  our  problem  this  way:  a volume  of 

oil  V per  unit  direction  perpendicular  to  x is  released 
*t  point  x - 0.  What  will  be  the  maximum  abscissa  x^ 
reached  by  the  oil  when  it  has  come  to  a standstill? 

We  will  Just  state  the  one-dimensional  results  since 
their  derivation  is  quite  similar  to  that  of  the  two-dimen- 


■ 


sional  study  in  the  previous  sub-sections.  A more  detailed 

analysis  is  given  in  Sub-Appendix  C-5.  The  first  statistical 

moment  <x  > of  the  oil  pool  final  size  is  found  to  be: 
max 


<x  > - 
max 


<h  (0> 


whereas  the  second  moment  a defined  as: 

x 


X - <X  5 

max  max 


appears  to  be: 


<x  >/ST 
~roax 


, r™4*  r<x  > nr  «u~ 

7 J L~  j L hJ 


dt  (31) 


Here  again,  ^ (x)  and  <h”>  are  related  to  t|»(x)  and  o. 

r 

through  the  relations: 


r>  - - \/5T 

V 


i|T (x)  - 


- 1|/) 

+ w 

o 


(u2  - j)  erf  (u) 


f 2 *o  - *\  . . *o  + ** 

\ ♦•+  V 


where:  - tjj(O) 

<J>  - <Kx) 
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3.6  General  Mathematical  Results 


In  this  section,  we  have  been  able  to  give  an  analyti- 
cal expression  for  the  first  two  statistical  moments  of  the 
oil  slick  final  size,  resulting  from  a spill  of  volume  V. 
Thus,  when  the  oil  spreading  is  either  one-dimensional 
(with  a privileged  direction)  or  two-dimensional  (no  priv- 
ileged direction) , the  assumptions  made  in  Section  2 enabled 
us  to  derive  the  ensemble  average  value  of  the  pool  size 
and  Its  standard  deviation,  as  well  as  a "dissymmetry  factor" 
M characterizing  the  aspect  of  the  slick  (In  the  two-dimen- 
sional case  only) . 

These  quantities  are  all  determined  in  terms  of  0^, 
the  variance  of  the  roughness  height  h,  and  the  auto-correla- 
tion function  of  h (either  H'(t)  or  4'(5,n). 


Namely: 


where  ^ (x)  is  obtained  by  equation  (15)  of  Sub-Appendix  C-3. 
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max 
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'1 
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<ls 

max 


<S  >]2> 
max 


>//T 


> 


2 -dim. 
case 


•f 


I 


k i 


i 


(35) 


where  ¥ is  obtained  by  Eq.  (18)  of  Sub-Appendix  C-3. 

However,  we  have  to  point  out  again  that  these  results 
do  not  constitute  the  unique  answer  to  the  general  problem 
of  oil  spreading  over  an  ice  surface  since  the  basic  assump- 
tions made  in  Sub-section  3.2  are  quite  limiting.  Never- 
theless, they  can  provide  us  with  a first  understanding  of 
the  possible  effects  of  an  oil  slick  in  an  Arctic  environment. 


4.  APPLICATION  - SAMPLE  CALCULATIONS 

i 

In  this  section,  the  general  mathematical  results  of  the 
last  section  will  be  applied  to  a certain  number  of  specific 
cases  of  ice  topography  where  the  corresponding  statistical  field 
data  is  available. 


i 

•I 


4 

| 


140 


Both  In  the  one  and  two-dimensional  cases  (for  privileged 
and  no  privileged  direction  of  spreading  respectively),  the 
statistics,  i.e.,  the  first  and  second  moments  of  the  final 
size  of  the  oil  pool,  have  been  determined  in  terms  of  o^,  the 
variance  of  the  roughness  and  tp,  its  auto-correlation  function. 
We  recall  here  that  h stands  for  the  variations  in  altitude  of 
the  ice  surface  with  mean  value  0. 


ah  - <h2>  (36) 

<Kt)  ■ <h(x)h(x+t)>  in  the  1-dim.  case  (37) 


'f(5»n)  • <h(x,y)h(x  + £,y  + n)>  in  the  2-dim.  case  (38) 


Then  within  a certain  region  of  the  Arctic  Slope  that  can  gen- 
erally be  regarded  as  large  compared  to  the  area  of  an  eventual 
oil  slick,  the  ice  topography  obeys  the  statistics  expressed 
either  by  Eqs.  (36)  and  (37)  or  by  Eqs.  (36)  and  (38)  when  the 
oil  spreading  phenomenon  is  one  and  two-dimensional  respectively. 

Now,  in  order  to  give  a numerical  answer  to  a specific  oil 
slick  problem,  we  have  to  know  and  either  ^ or  f character- 
izing that  part  of  the  Arctic  Slope  where  the  spill  occurs,  and 
that  are  the  results  of  field  observations.  However,  the  statist- 
ical properties  of  the  permafrost  relief  are  generally  given  in 
the  form  of  a power  density  spectrum  of  the  ice  surface  (Ref.  [3]). 
We  must  therefore  use  that  kind  of  field  data  Instead  of  the 
auto-correlation  function  itself. 
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where  U is  a frequency  vector  of  components  and  and  X any 
vector  of  the  plane  xy  of  components  £ and  n« 


Therefore,  It  la  clear  that  the  information  we  need  about 
the  Ice  surface  statistics,  i.e.,  c and  the  auto-correlation 
function  of  h,  Is  contained  In  the  power  density  spectrum  P. 

The  problem  Is  now  to  build  a continuous  function  P fitting  the 
observed  power  spectrum  as  well  as  possible.  This  is  the  object 
of  the  next  sub-section. 


1 


I 


1. 1 


i 1 


: 


I 


4.1  Modelling  of  the  Power  Spectra 

In  order  to  make  the  mathematical  analysis  of  Section  2 
applicable  to  any  part  of  the  ice  in  the  Arctic  Ocean, 

we  have  to  determine  the  parameters  necessary  to  describe 
the  most  important  features  of  the  ice  surface  statistics. 

A general  answer  to  this  modelling  problem  cannot  be  given 
here,  even  by  increasing  the  number  of  such  parameters 
since  the  actual  data  is  far  from  being  complete  and  re- 
liable. Moreover,  another  problem  is  raised  by  the  two- 
dimensional  case  of  oil  spreading  since  no  corresponding 
field  data  - two-dimensional  functions  Y(£,n)  - is  available 
to  date.  Therefore,  in  order  to  obtain  numerical  orders  of 
magnitude,  we  have  assumed  that  a two-dimensional  auto- 
correlation function  V(£,n)  can  be  regarded  as  the  product 
of  two  one-dimensional  auto-correlation  functions.  Namely: 


Y«.n> 


0.  <M5)  Mb) 


(41) 
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There  is  no  evidence  at  all  that  this  is  correct  and  It 
only  appears  as  a fair  and  indispensable  assumption  to  treat 
the  two-dimensional  case.  Therefore,  since  the  latter  can 
be  derived  from  one-dimensional  field  data,  this  sub-section 
will  only  be  concerned  with  the  modelling  of  the  one-dimen- 
slonal  power  spectra  that  directly  come  out  by  Fourier 
analyzing  laser  profilometer  traces  such  as  those  mentioned 
in  the  Introduction.  The  two-dimensional  model  will  only 
be  stated  following  the  above  assumption. 

There  is  no  doubt  that  modelling  P(f)  or  <Kt)  Is 
equivalent  since  the  power  density  spectrum  and  the  auto- 
correlation function  are  uniquely  related  through  a Fourier 
transformation.  In  fact,  from  the  few  data  examined  - 
both  ijs  and  P -,  it  turns  out  that  the  shape  of  the  spectra 
is  quite  varying  thereby  making  the  choice  of  parameters 
difficult,  whereas  the  form  of  ^(T)  seems  to  remain  pretty 
much  the  same;  it  is,  therefore,  easier  and  perhaps  more 
reliable  to  propose  a model  for  <Kt)  and  then  take  its 
Fourier  transform.  Since  this  function  presents  observed 
damped  oscillations  around  a zero  value  and  starts  at 
ij»(0)  - 0^  (by  definition  of  i|<),  we  think  that  two  parameters 
k and  T are  sufficient  to  set  up  a consistent  model.  Namely: 

T 

*(t)  - cos  2tt  £ (42) 

1 + kt  1 


where  k is  the  Inverse  of  a length  squared  and  T has  the 
dimension  of  a length. 
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However,  as  we  mentioned  at  the  beginning  of  this 
section,  this  is  not  the  form  under  which  field  measure- 
ments generally  appear.  Thus,  we  have  to  derive  a model 
for  the  power  spectrum  P(f).  This  can  be  done  by  taking 
the  inverse  Fourier  transform  of  the  model  for  i|>  and 
evaluating  k and  T from  the  power  spectra.  It  follows 
that: 


P(f)  . — exp  (-  — ^coshf  — f > 1/T 
A \AS  \T  AJ 


P(f)  “ 


TO. 


A 


exp  ( fe)  coshCf) f < 1/1 


> 


P(f)  - p(-f) 


(A3) 


There  exists  a symmetry  with  respect  to  the  origin  and  it 
is  therefore  sufficient  to  consider  only  the  positive  fre- 
quencies for  which  a peak  in  the  power  spectrum  appears  at 
F - 1/T.  In  that  case,  the  normalization  of  the  spectrum 
is  such  that: 


P(f)  df 


(AA) 


which  is  equivalent  to  multiplying  P(F)  by  a factor  two 
in  relation  (A3) . 


1A5 


This  relation  constitutes  the  two-parameters  raodel^ 
we  propose  for  the  one-dimensional  power  spectra.  When 
considering  specific  ice  surface  statistics,  the  parameters 
k and  T can  be  determined  by  having  the  model  fit  exactly 
the  observed  intensity  and  frequency  position  of  the  peak 
on  the  power  spectrum.  In  other  words,  if  f is  observed 
position  of  the  peak,  we  have: 

T - l/f  (45) 

no 

and  k is  the  solution  of  the  following  transcendental 
equation: 

p(f  ) ■ exp  (-  ^L\coahf2l_\  (46) 

no  A \ T A)  \TAS 

Figures  4 to  6 show  one-dimensional  spectra  modelled  in 
this  way.  Each  figure  is  relevant  to  different  ice  surface 
statistics,  i.e.,  different  k's  and  T's;  the  solid  line 
shows  the  proposed  model  whereas  the  black  squares  represent 
the  corresponding  field  measured  spectrum.  It  appears  that 
the  decaying  of  the  spectrum  is  acceptably  verified  but, 
unfortunately,  the  lower  frequency  region,  beyond  the  peak 
is  not  very  well  matched.  As  the  peaks  of  the  power  density 
spectra  play  a central  role  in  our  analysis,  a point  per- 
taining to  their  exact  definition  has  to  be  raised  here. 


*Oh  is  not  a parameter  in  that  model  since  J P(F)  df  ■ 
for  any  k's  and  T's  0 
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Aa  a matter  of  fact,  the  behavior  of  the  power  spectra 
for  very  low  frequencies  (where  the  peaks  lie  most  often) 
remain  unclear.  These  data  generally  come  out  by  Fourier 
analyzing  field  measurements  of  the  ice  surface  variations 
in  altitude.  We  must  notice  that  the  recording  device 
does  not  have  an  Infinite  analyzing  band-width  and  that  in 
order  to  remove  the  variations  in  altitude  of  the  plane 
which  carries  the  recorder  the  data  signal  is  filtered 
through  a hagh-band  pass  filter.  Therefore,  we  cannot 
rely  on  the  information  at  low  frequencies  provided  by 
such  methods,  and  it  would  be  necessary  to  set  up  a physical 
argument  to  account  for  the  general  aspect  of  the  power 
spectra  in  the  low  frequencies.  This  requires  a thorough 
understanding  of  the  ice  formation  and  life  and  has  not 
been  done  to  date.  Our  modelling  is  therefore  restricted 
to  such  field  measured  power  spectra  that  exhibit  a well 
defined  single  predominant  peak,  almost  regardless  of  their 
behavior  for  frequencies  less  than  the  peak  frequency. 

Let  us  now  state  the  resulting  model  for  a two-dimen- 
sional statistics.  Following  the  assumption  we  made  earlier, 
we  write: 

v«,n)  - oh  *2(o)  (47) 

where  ^ and  ^ are  one-dimensional  auto-correlation  func- 
tions, hence  modelled  by  relation  (42)  with  parameters  k^, 
and  k2  respectively.  It  follows  that  the  corresponding 
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two-dimensional  power  spectrum  is  P(U)  such  that: 


P(U)  - ah  P1(u1)  p2(u2) 


where  and  U 2 are  the  components  of  the  frequency  vector 


U and  P^  and  P2  are  the  two  one-dimensional  models  cor- 


responding to  the  parameters  k^,T^  80 ^ ^2’^2  resPect^ve^y' 
Namely: 


W 


P2(U2) 


P1(U1) 


p2(u2) 


/ 2uU, 


» / " '1\  , / 2tt 

exp  [ 1 cosh 


^1  V*1 


T /k 

V 1‘ 


_ / 2ttU 

exp  ( | cosh 

^1  \^2 


(T2^) 


if  UL,U2  > 1/T 


i_  exp  / ^\C0,,/!!!A 

vAi  \*J 


exp  - 


if  UltU2  < 1/T 


P(-U)  *=  P(U) 


Here  again,  P(£)  is  symmetric  with  respect  to  tne 
origin  in  the  frequency  plane.  Hence,  one  needs  only 


consider  a half-plane  (for  instance  > 0),  normalizing 
the  power  density  spectrum  so  that: 


(50) 


which  is  equivalent  to  multiplying  P(U)  as  expressed  in 
Eqs.  (48)  and  (49)  by  a factor  2. 

The  overall  aspect  of  such  a power  spectrum  is  given 
in  Fig.  7.  Two  peaks  appear,  corresponding  to  the  two 
peaks  of  P2(U2>.  As  was  said  earlier,  no  comparison  with 
field  data  is  possible  due  to  the  lack  of  the  latter. 

In  this  section,  a simple  two-parameters  model  has 
been  proposed  for  the  one-dimensional  power  spectra  and 
is  in  good  agreement  with  observations  in  so  far  as  a real 
predominant  peak  appears  in  the  measured  spectra.  The  two- 
dimensional  model  has  been  purely  hypothesized  as  being  a 
cross  product  of  two  one-dimensional  models.  The  next  sub- 
section will  use  these  models  to  develop  a concept  of  "most 
probable  pocket"  in  the  ice  relief. 


Concept  of  Most  Probable  Pocket 


In  this  sub-section,  both  in  the  one  and  two-dimensional 
cases,  a "most  probable  pocket"  will  be  defined  as  the  "most 
apparent  periodical  relief"  in  the  ice  surface;  mathemati- 
cally, it  will  be  identified  with  the  negative  part  of  the 
main  tern  in  the  Fourier  series  expansion  of  the  ice  surface. 
Moreover,  being  the  most  characteristic  feature  of  an  ice 
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surface,  this  most  probable  pocket  will  be  used  to  non- 
dimensionalize  our  further  results. 

As  soon  as  we  have  the  power  density  spectrum  P(w)  of 
the  function  h,  from  harmonic  analysis  we  know  that: 


h(x) 


F(w)  cos  (wx)  dw 


with  2P(w)  dw  - | F(w) | ‘ 

w ■ 2 iff 


(51) 


It  is  then  clear  that  we  can  approximate  this  function 
by  an  infinite  Fourier  cosine  series  which  is  nothing  more 
than  replacing  the  surface  contained  between  the  power 
spectrum  curve  and  the  frequency  axis  by  an  infinite  num- 
ber of  rectangles  of  width  Aw^  and  height  P(wn)  where  wq 
is  a circular  frequency  and  P(wr)  the  value  of  the  power 
spectrum  at  that  frequency. 

Thus: 


00 

h(x)  - E /2P(w  ) Aw  cos  (w  x - i|/)  (52) 

q n n n n 

in  the  one-dimensional  case 

by  the  same  token: 

00 

h(X)  = E /ZP(U)  U"  cos  (U  .X  - to)  (53) 

— -n  “ti  — rn 

in  the  two-dimensional  case 

T 

See  Ref.  4. 
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In  Chcse  equations,  the  "approximately  equal"  symbol 
tends  to  an  "equal"  symbol  when  Aw^  - AU^  in  the  two-dimen- 
sional case  - tends  to  zero,  that  is  to  say  when  the  sum  is 
replaced  by  an  Integral. 

The  amount  of  randomness  characteristic  of  any  given 
realization  of  the  ice  surface  among  a particular  statistics 
appears  in  the  set  of  phases  that  are  regarded  as  com- 
pletely random  variables. 

Now,  in  so  far  as  the  power  spectrum  decays  fast  enough 

when  the  frequency  Increases  and  decreases  from  its  value 

at  the  peak,  there  is  no  doubt  that  the  term  corresponding 

1 2 + 

to  the  frequency  position  of  the  peak,  w - U , IT  in 

no  —no  —no 

the  two-dimensional  case  - in  the  Fourier  series  is  pre- 
vailing over  the  other  ones. 

Figures  4 to  6 show  that  such  a fast  decaying  is  con- 
firmed by  field  data  when  w increases,  but  this  is  not  the 
case  when  w decreases  since,  as  was  already  mentioned,  this 
low  frequency  part  of  the  spectra  is  not  quite  reliable. 
However,  it  seems  intuitively  reasonable  to  expect  that  the 
zero  frequency  values  of  the  spectra  have  to  be  zero;  in 
other  words,  in  the  Fourier  expansion  of  h,  terms  with  an 
infinite  wave  length  should  have  zero  amplitude  because 
otherwise,  this  would  mean  that  there  exists  points  in- 
finitely far  from  the  origin  of  the  coordinates  with  non- 
zero roughness  height.  This  is  intuitively  non-acceptable. 

A rigorous  argument  can  only  come  out  from  a detailed 
discussion  of  the  ice  relief  formation  and  evolution.  This 
is  beyond  the  scope  of  this  analysis  and  we  must  be  content 


t 1 2 

U , U“  refer  here  to  the  two  peaks  of  a two-dimensional 
-no  -no  r 

spectrum. 
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with  the  assumption  that  a real  peak  exists  in  every  spectra 
and  its  frequency  position  and  intensity  explicitly  appear 
on  the  field  data. 


Thus,  it  follows  that  an  even  cruder  approximation  for 
h would  be: 


h(x)  a /ZoT  cob  (w  x - ) (54) 

n no  no 

in  the  one-dimensional  case 


and 


h(X)  s /2^ 


cos  (U 

L. 


-no 


> 


+ cos  (U  .X  - 
-no  — 


(55) 


in  the  two-dimensional  case 


For  the  one-dimensional  case,  we  can  look  at  Eq.  (52)  as 

expressing  h in  terms  of  a superposition  of  sinusoids  of 
2tt 

period  — and  intensity  /2P(w  ) Aw  . Each  sinusoid  will 
n 

contribute  - in  proportion  to  its  amplitude  - to  form  the 
actual  aspect  of  a typical  one-dimensional  profilometer 
trace  verifying  the  statistics  expressed  in  the  power  spec- 
trum P.  Then,  Eq.  (54)  tells  us  that  the  greatest  ampli- 
tude sinusoid,  and  therefore  that  which  most  influences 
the  shape  of  the  real  pockets,  corresponds  to  the  frequency 


of  the  peak.  Hence  the  notation  "most  probable  pocket" 
refers  to  this  particular  term  of  the  Fourier  expansion 
of  h.  Tills  privileged  pattern  in  the  ice  relief  has  a 
length  Lmp  equal  to  one  period  T of  the  sinusoid  and  a 
"volume"  per  unit  width  Imp  equal  to  the  area  contained 
between  a negative  arch  of  the  curve  and  the  axis,  namely: 


Imp 


2tt 

w 

no 


T 


(56) 


Imp 


/2a,  cos  (w  x)  dx 
n no 


negative  integrand 


Or,  integrating: 


(57) 


Emp 


/i3h 


T /to” 
n 

TT 


(58) 


By  the  same  token,  for  the  two-dimensional  case, 

Eq,  (53)  is  an  expansion  of  the  ice  surface  in  terms  of 
sinusoidal  surfaces,  the  most  apparent  pattern  corresponding 
to  the  two  peaks  of  frequency  vectors  ■ (UnQ^,  ^no2^ 
and  Jj£o  - (Unol,  - U^)  as  in  Eq.  (55),  (see  Fig.  7). 
Analogous  to  the  one-dimensional  case,  the  "most  probable 
pocket"  has  here  a surface  S and  a volume  corresponding 
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to  on  increase  of  (Uro^)  * in  x and  (U^)  * in  y* 
Or  more  precisely: 


1 period  in 
1 period  in 
negative  Integrand 


Upon  Integration,  this  yields: 


# 


Vmp  - 


“nol  Uno2 


(61) 


Since  they  express  the  most  Important  feature  of  the 
power  density  spectra  - namely,  their  peak  -,  the  "most 
probable  pocket"  dimensions  will  provide  a good  reference 
to  non-diraensionalize  our  further  results.  This  concept 
of  most  probable  pocket  has  also  the  advantage  of  giving, 
in  a simple  way  (through  Eqs.  (52)  and  (53),  a rough  esti- 
mate of  what  the  ice  surface  looks  like. 
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4.3  General  Mathematical  Results  of  Section  3 
Using  the  Modelled  Spectra  of  4.1 


With  Eqs.  (34)  In  the  one-dimensional  case  (35)  In 
the  two-dimensional  case,  we  have  obtained  a relation 
between  the  first  and  second  statistical  moments  of  the 
oil  pool  final  size,  the  auto-correlation  function  of  h 
and  the  parameter  of  the  ice  surface  statistics.  The 
purpose  of  this  section  is  to  derive  more  simple  relations, 
using  modelled  spectra,  and  non-dimens ionalize  these  re- 
sults with  rtspect  to  the  "most  probable  pocket".  Thus, 
it  will  be  shown  that  in  their  dimensionless  form,  the 
second  moment  is  proportional  to  the  first  one  through  a 
factor  8 only  dependent  on  the  non-dimensional  group 
kT2  (or  k^k^  *n  t*ie  two-dimensional  case)  for  suit- 

able values  of  this  group. 

4.3.1  One-Dimensional  Case 

We  recall  here  that  if  <x  > is  the  final 

max 

dimension  of  the  oil  pool  and  a ^ its  variance,  we 
have,  from  the  mathematical  analysis  of  Section  3: 


i|/”(T)  being  obtainable  from  <KT)  defined  as: 


Wt)  ■ — 7 cos  2ir  £ 

1 + kT4 


(64) 


Through  the  relation  (see  Sub-Appendix  C-4) . 


f (T)  - ^(W 


2lr(*02 


ab 


(65) 


We  can  notice  that,  for  k’s  and  T's  such  that 
approaches  zero  within  one  period  T - 
or  one  most  probable  pocket  length  Imp  - Eq.  (63) 


We  can  r 

L°»  2!J 


becomes: 


(66) 


Since  Che  integrand  in  relation  (66)  is  a con- 
tinuous function  of  T , we  can  define  T such  that: 

o 


0^  **  8tt  Lrap 


r<x  > 

max  _ „ 

1 

H 

0 

1 

1 *”* 

! 1 

1 

isSj 

1 

Li 

/h  27j 

(67) 


Mow,  from  Sub-Appendix  C-4,  it  it  clear  that 
will  be  positive  between  T * 0 and 


£m.  iT 

Lah  2"J 


T T T 

T ■ and  negative  between  t - -j-  and  x ■ j;  hence: 


0 < T < T/4 
o 


(68) 


In  so  far  as  <x  > is  large  as  compared  to  -f, 
max  4 

we  can  rewrite  Eq.  (67): 


- 8ir  . . 

0 ■ — Lmp  a <x  > 

x /2  maX 


(69) 


with 


a ■ 


♦‘(y  1. 

°h  2* 


(70) 
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Taking  the  square  root  of  Eq.  (69),  it  appears 

that  the  standard  deviation  of  the  oil  pool  final 

1/2 

size  is  proportional  to  <x  > . This  is  in  ac- 

max 

cordance  with  the  following  physical  argument.  Let 
V be  the  volume  of  a pocket  per  unit  width,  and 
let  us  split  it  into  a "most  probable  pocket"  volume 
plus  a fluctuating  quantity  V': 

AV  ■ Lmp  <h  > + V'  (71) 


If  is  the  variance  of  AV  and  Oy,  that  of 
V',  we  have  obviously: 


AV 


(72) 


Now  that  the  total  volume  V of  a spill  is  such 

that: 


V - N AV 


(73) 


Furthermore,  all  these  pockets  are  statistically 
independent  if  we  assume  that  the  correlation  for  h 
goes  to  zero  within  an  Interval  Lmp. 


T 


Lmp  is  defined  in  (56). 
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Therefore 


N o 


AV 


N a 


r 


(74) 


where  is  the  variance  of  the  volume  occupied  by 
N pockets.  Since  <x  > is  always  of  the  form 


<x  > 
max 


■ S Bp 


(75) 


we  have  from  (74): 


O 


V 


<x  > 


max 


(76) 


1/2 

1/2 

is  then  proportional  to  <x  > 

max 

as  expressed  in  Eq.  (69). 

Now,  let  us  show  that  a defined  in  relation 

(70)  is  independent  of  <xmax>  for  sufficiently  large 

values  of  <x  >. 

max 
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Equation  (66)  can  be  rewritten: 


i<x  >/^r  - t] 


max 


[<x  >//!  - T ] 

max  o 


(^'(T)/oh  - ij-)  dT  (77) 


the  ratio 


<X  >1/1  - T 
max 

<x  >1/2  - t 
max  o 


can  be  regarded  as  constant 


in  so  far  T - T remains  small  as  compared  to  <x  > 
o max 

when  t varies  from  0 to  Imp.  In  other  words,  since 

T is  of  order  mp/4,  a independent  of  <x  > is  a 
o max 

reasonable  approximation  when  the  number  of  "most 
probable  pockets"  involved  in  the  oil  spill  is  large 
(of  order  10  or  greater). 

2 

We  can  also  shew  that,  when  kT  exceeds  a number 
to  be  determined,  a only  depends  on  the  non-dimensional 
group  kT  : 


iT(t) 


-°F? 


(*)  - 


2tt(*o2  - <J>2) 


1/2 


abexp 


(a2  + b2)i|»  - 2ab^ 
20|>o2  - *2) 


dadb  (78) 


'if 
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Furthermore,  it  appears  that,  setting 


Non-dimensionalizing  Eq.  (69)  with  respect  to 
the  "most  probable  pocket"  length  Lmp,  we  obtain: 


with 


0 - 


(82) 


(83) 


We  now  determine  the  universal  curve  6 versus  kT 
for  k and  T's  that  verify  the  assumptions  we  have 
made  in  the  previous  derivations. 


The  only  one  concerned  with  k and  T tells  us 


must  decay  to  zero  within  a period 


t 


T ■ Lmp.  From  Sub -Appendix  C-4,  we  know  that  this 
is  equivalent  to  ip(x)  going  to  zero  within  a period 
T;  if  we  assume  that  i|/(T)  ■ --^p-  *8  a good  approxi- 
mation for  <p(T)  equals  zero,  it  follows  that: 


iK0)/<KT)  >,  50 


(84) 


We  could  have  taken  the  same  condition  within  two  or 
more  periods  T;  the  dependence  of  a on  kT^  only  would 
not  have  been  changed,  but  a independent  of  <x  > 
could  be  questioned.  However,  the  decaying  to  zero 
within  a period  seems  to  be  in  good  accordance  with 
field  data. 
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■ 


or: 


kT2  * 50  (85) 


2 

We  will  therefore  compute  the  curve  0 (kT  ) 

2 

for  kT  2 50.  This  is  represented  on  Fig.  11. 

2 

We  can  notice  that,  when  kT  goes  to  infinity, 

6 goes  to  zero.  This  is  not  surprising  since, 
when  T ■ ®,  a equals  a constant  divided  by 
infinity  and  when  k ■ ®°,  a equals  zero  times  a 
constant;  in  both  cases,  6 is  therefore  zero. 

Knowing  the  k and  T of  a given  ice  surface 

statistics,  we  are  now  in  a position  to  obtain 

the  corresponding  0 (on  Fig.  11)  and,  with  Eq.  (82) 

have  got  a single  relationship  between  the  two 

first  statistical  moments  of  the  oil  pool  size 

x , in  their  non-dimensional  form, 
max 
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4.3.2  Two-Dimensional  Case 

We  recall  that  the  two-dimensional  size  of  the 

oil  pool  is  S and  its  variance  a such  that: 
max  s 


<S  > - ^ 

max  /ST 


(86) 


O - 32  TT 

8 


<x  >//2  <y  >/S2 
max  ymax  p<x  > 

/ / BH 


<y  > 
7 max 

^"(C.n)  l 

- n 

_/2 

. 0h  2’_ 

dCdn  (87) 


where  V~($,r\)  i8  obtainable  from  4'(£,n)  defined  as: 


ne.n)  » 2 


(l  + kjOU  + k2n  ) 


cos  2tt^—  cos  2tt~— 
*1  X2 


(88) 
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Going  through  the  same  kind  of  derivations  as 

ve  did  for  the  one-dimensional  case,  and  under  the 

same  conditions  for  k's,  T's  and  S , we  will  ob- 

max 

tain  comparable  results. 

Thus,  for  k^,  k2>  and  T2  such  that  ¥ (£,n) 

goes  to  0 within  the  area  X T2  and  sufficiently 

large  <S  >,  we  can  write 
max 


o1'2 

8 

Smp 


(89) 


2 2 

0 being  only  a function  of  k k.  T.  T«  (if  we 

2 2 1 *•  L L 22 
neglect  kj^  + K2  T2  compared  to  k^  k2  T2  ). 

0 versus  dimensionless  group  has  also  been  deter- 
mined on  Fig.  13,  and  one  notices  that  for  ¥ (£,n) 

2 2 

to  go  to  zero  within  the  area  X T2,  k^  k2  T2 
has  to  be  greater  than  2500  if  we  keep  the  same  stan- 
dards as  in  the  one-dimensional  case  for  the  meaning 
of  "goes  to  zero". 

Therefore,  knowing  k^,  k2>  and  T2  of  the  two- 
dimensional  statistics  we  are  considering,  we  are  in 
a position  to  obtain  the  corresponding  0 and  thereby, 
we  have  got  a simple  y power  law  relationship  be- 
tween the  non-dimensional  first  and  second  statistical 
moments  of  the  oil  pool  size. 
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4.  4 Results 


Having  three  experimental  power  density  spectra  coming 
from  different  parts  of  the  Arctic  Ocean,  we  have  attempted 
to  insert  them  in  our  parametered  model  described  In  Section 
4.1.  We  obtained  the  following  values  for  the  parameters 
k and  T (see  the  modelled  spectra  on  Figs.  4 to  6) . 


T 

k 

°h 

115  ft. 

0.00378  ft.“2 

2.295  ft.2 

200  ft. 

0.002786  ft.-2 

1.016  ft.2 

540  ft. 

0.04075  ft.-2 

0.2033  ft.2 

© 

© 

© 


(90) 


Then,  on  Fig.  10,  plotting  the  relationship  between 
the  two  first  moments  of  the  oil  pool  size  in  the  one- 
dimensional  case  for  any  of  the  three  above  statistics, 
we  have  compared  the  simple  ^ power  laws  given  by  Eq.  (82) 
with  the  corresponding  0's  (solid  lines)  to  points  obtained 
from  the  general  mathematical  model  - Eq.  (34)  - for  dif- 
ferent volumes  released  V.  The  agreement  is  extremely 
good  for  any  statistics,  as  well  as  it  is  on  Fig.  12  for 
the  two-dimensional  case;  in  order  to  build  two-dimensional 
statistics,  we  have  made  the  following  combinations  of  one- 
dimensional  statistics:  ® + ©>©  + ©»©  + CD 
® + © of  Table  (90). 

+ - 1 n ^ 

We  notice  that  kT  is  always  larger  than  50  for  all  three 
statistics  which  allows  us  to  use  Section  4.3  analysis. 
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Therefore,  In  so  far  the  conditions  over  k and  T is 

2 1 
verified  (kT  ^ 50  for  instance),  this  j power  law  re- 
lationship between  the  oil  pool  final  size  and  its  stan- 
dard deviation  appears  to  be  a very  good  approximation  to 
solve  our  problem,  both  in  the  one  and  two-dimensional 
case. 

It  might  be  interesting  to  f.how  a rough  estimate  of 
what  these  numerical  results  can  be.  Let  us  take  a volume 

3 

V * 1000  ft.  per  foot  width  in  our  one-dimensional  case; 
this  yields  the  following  figures  for  each  of  the  statistics 
of  Table  (90): 


<x  > 


max 


1458  ft. 


a 1/2  - 304  ft. 

x 


<x  > 
max 


2487  ft. 


480  ft. 


<x  > 
max 


5560  ft. 


675  ft. 


As  to  the  two  dimensional  case,  let  us  take  a volume 
5 3 

V ■ 5 X 10  ft.  which  appears  to  be  a reasonable  size, 
even  in  case  of  a supertanker  crash,  then,  we  obtain  for 
the  following  combinations  of  one-dimensional  statistics: 
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ip?*'-  v 


'ii 


Such  results  for  0 and  o provide  good  evidence  that 

X 8 

our  assumptions  of  Eqs.  (11)  and  (7,  in  Sub-Appendix  C5) 

are  valid.  We  recall  that  in  these  assumptions,  the 

- 1/2  1/2 

average  of  h over  a length  a or  area  a was  identi- 

■ ■ m — — ■ i X S 

fled  to  h (O , the  average  of  h over  the  whole  reallza- 

1/2  1/2 

tlon  of  the  surface,  on  the  ground  that  or  Og  was 

large  enough. 

Since  the  correlation  for  h approaches  zero  within  one 
"most  probable  pocket"  dimension,  a length  I-mp  - or  an  area 
Smp  in  the  two-dimension  case  - of  the  ice  surface  realiza- 
tion can  be  considered  as  containing  all*  the  statistical 
characteristics  of  the  entire  surface.  Therefore,  the 
average  value  of  h over  a length  Lmp  - or  area  Smp  should 
not  differ  very  much  from  <h  (£)>. 

Hence,  it  comes  out  that  for  V such  that  the  standard 

deviation  of  the  oil  pool  is  much  less  than  a "most  probable 

pocket",  the  assumptions  are  no  longer  verified,  thereby 

making  the  lower  part  of  the  curves  0 versus  <x  > - or 

x wax 

0 versus  <S  > - less  trustworthy, 
s max  1 


5.  CONCLUSION  - SUMMARY  - DISCUSSION 

In  this  study,  the  two  first  statistical  moments  of  the 
oil  pool  resulting  from  a spill  of  volume  V have  been  deter- 
mined in  terms  of  O^,  the  variance  of  the  roughness  height  h 
and  the  auto-correlation  function  of  h,  both  when  the  spreading 
is  one-dimensional  (one  privileged  direction)  and  when  it  is 
two-dimensional  (no  privileged  direction). 
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These  results  have  then  been  developed  and  simplified  in 
compliance  with  a two-parameter  (k,T)  model  proposed  for  one- 
dimensional statistics  as  long  as  their  power  spectrum  exhibits 
a single  predominant  peak.  If  P(f)  is  the  modelled  spectrum, 
we  have: 


^exp  (-2StW2H-N 

A \A/  \tAJ 


if  f > l/T 


^exp  /.  2i_'\c<)Sh(' m\ 
A \iAj  \A  J 


if  f < l/T 


T is  nothing  more  than  the  inverse  of  the  frequency  position  of 
the  peak  and  k is  determined  by  the  intensity  of  the  spectrum 
at  that  peak,  through  the  transcendental  equation: 


p,lv  2TOh  7 2ir  \ 

P(  ) . exp cosh 

T A \ TAJ 


Two-dimensional  statistics  have  been  regarded  as  character- 
ized by  a product  of  two  one-dimensional  power  density  spectra, 
therefore,  their  modelling  is  quite  similar  to  that  of  the  one- 
dimensional  case. 
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K "most  probable  pocket"  that  could  be  physically  defined 
.a  the  'Wet  apparent  periodical  relief”  In  the  1c.  surface  has 
been  Introduced  both  In  the  one  end  two-dimensional  case.  Math- 
ematically, It  Is  Identified  to  the  main  term  In  the  Fourier 
series  expansion  of  the  ice  surface. 

If  ± is  the  frequency  position  of  the  peak  in  a one-dimen- 

T j /tt  ,t  \ /’ll  - U „)  the  positions 

sional  spectrum  and  (Unol»  un02;*  1 nol*  no2 

of  the  peaks  n a two-dimensional  spectrum  (see  Fig.  7),  the 
dlmenelo-  f this  "most  probable  pocket"  are  Its  size  and 

volume : 


Length  Lmp  “ T 

"volume"  £mp  ” « 


area  Smp 


volume  Vmp 


T1T2 


4tt 


Unol  Uno2 


8/2 o7 


U , U , 
nol  no2 


two  dim.  case 

(95) 


The  first  statistical  moment  of  the  final  size  of  the  oil 
pool  - namely  x^  in  the  one-dimensional  case  or  in  the 

two-dimensional  one  - has  been  determined  and  then  non-dimen- 
sionalised  with  respect  to  the  "most  probable  pocket"  size: 
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(96) 


<x  > 
max 

Lrap 


<S  > 
max 


Smp 


V 

Lmp 


V 

Smp 


\/lF-  JL. 

V ®p  /jf 


V 

Vmp 


ir 


4 

371 


A simple  power  law  relationship  between  the  non-dimen- 
sional first  and  second  moments  of  the  pool  size  has  been  de- 
veloped with  a slope  6 - in  the  two-dimensional  case  - 

2 2 2 

depending  only  on  the  dimensionless  group  kT  - k^  k^ 

in  the  two-dimensional  case  - that  characterized  the  modelled 


spectra. 

of  S : 
max 


Namely,  if  o is  the  variance  of  x and  a that 

•'  9 v mov  a 


max 


1/2 


x 

L mp 


■ 0 


>Y 12 


see  Fig.  (10)  and  (11) 


(97) 


1/2 

as 

Smp 


>\l/2 


see  Figs.  (12)  and  (13) 


(98) 


Therefore,  suppose  we  are  given  the  volume  V of  oil  released 
and  the  parameters  which  characterize  the  region  of  the  Arctic 
Ocean  where  the  spill  takes  place:  k and  T in  the  one-dimen- 
sional case  or  k^,  k^,  and  T2  in  the  two-dimensional  one, 

along  with  the  variance  of  the  roughness  height;  we  are  then 
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la  a position  to  predict  the  final  size  of  the  oil  pool  as  well 
as  Its  standard  deviation.  Furthermore,  In  the  two-dimensional 
case,  an  Idea  about  the  dissymetry  of  the  pool  can  be  gotten 
from  the  quantity  y,  the  ratio  of  the  dimensions  of  the  rectangle 
limiting  the  pool.  With  relation  (4) : 


U " 


<2ff/T2r  + 2 k2 


+ 2 kx 


1/2 


(99) 


We  have  to  notice  here  that  this  rectangular  shape 

x X y surely  does  not  represent  the  real  aspect  of  an  oil 
ni8X  max 

pool.  The  latter  is  much  closer  to  an  ellipse  - a circle  in  the 
case  of  a truly  isotropic  surface  - and  it  follows  that  the  size 
computed  under  this  assumption  can  result  in  a 13%  underestima- 
tion of  <x  > and  <y  >.  The  computations  relative  to  an 
ellipse,  although  quite  intricate,  can  be  performed  but,  in 
view  of  the  precision  obtained  on  the  statistical  field  measure- 
ments, we  did  not  deem  it  worthwhile. 

The  limitations  for  the  use  of  these  power  law  relationships 

between  the  first  two  moments  of  the  pool  size  (Eqs.  (97)  and 

2 

(98))  are  concerned  with  the  value  of  kT  and  <x  > - <S  > 
22  max  max 

and  k^  k2  T^  T2  in  the  two-dimensional  case  - namely: 


kT2  » 50  or  kx  k£  Tj2  T*  >,  2500 


(100) 
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<x 


<S  > 

max  / _max_  ^ 1Q 


Lmp 


Lnp 


(101) 


These  conditions,  described  in  some  detail  in  Sub-section 
4.3  do  not  influence  more  general  results  as  those  of  Eqs.  (34) 
and  (35)  which  are  valid  for  any  power  spectra  model. 

Also,  an  important  point  to  understand  is  the  actual  be- 
havior of  the  power  spectra  at  very  low  frequencies  since  it 
has  been  here  purely  hypothesized,  as  was  also  the  two-dimensional 
statistics,  for  there  is  no  evidence  at  all  that  a two-dimensional 
power  spectrum  should  be  the  product  of  two  one-dimensional  ones. 

To  conclude  this  study,  we  can  say  that  the  previous  analysis 
yields  simple  and  useful  results  as  to  the  prediction  of  oil 
spreading  over  a rough  surface  of  ice;  but,  since  the  modelling 
of  the  power  density  spectra  has  been  developed  from  a limited 
quantity  of  data,  a strong  experimental  ground  is  required  to 
completely  justify  this  proposed  model. 
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FIG.  I GENERAL  ASPECT  OF  ONE  REALIZATION  OF 
THE  PROFILE 


FIG.  2 AVERAGE  LONGITUDINAL  DIMENSION  OF 
A POCKET  ' 
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VERIFICATION  FOR  h NORMALLY  DISTRIBUTED 


FIG.  7 EXAMPLE  OF  TWO-DIMENSIONAL  POWER  SPECTRUM 


FIG.  8 ASPECT  OF  THE  AUTO-CORRELATION  FUNCTION  O 
h“  AND  h+ 
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THIS  AREA  EQUALS 
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G.  9 POSITION  OF  r0  WITH  RESPECT  TO  T 
(SECTION  m-  3 OF  THE  TEXT) 


FIG.  10  RELATIONSHIP  BETWEEN  <XMAX>  AND  o^,/2. 
ONE-DIMENSIONAL  CASE 


TWO-DIMENSIONAL  CASE 


F I* 


t | 


l 

\ 

t 
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SUB-APPENDIX  C-l 

In  this  appendix,  we  will  determine  the  "dissymmetry  factor" 
defined  as  the  ratio 


Xx 

11  - Tv 


(l) 


where  Xx  and  Xy  are  the  average  dimensions  of  a pocket  In  the  x and 
y directions  respectively. 

Let  us  consider  an  actual  laser  profilometer  trace  parallel  to 

the  x direction.  Let  r (x)  denote  the  elevation  of  ice  above  the 

x 

mean  level  < h(x)  > ■ 0.  Considering  an  ensemble  average  over  all 
the  traces  parallel  to  x,  we  can  write: 


< r 00  > - 0. 
x h 


(2) 


And: 


< rx(x  + O rx(x)  > 


*«.<>> 


(3) 


Using  our  third  hypothesis  concerned  with  the  way  pockets  are  filled 
by  oil,  we  can  see  that  Xx  corresponds  to  the  average  distance 
separating  two  zeros  of  the  function  rx(x)  (see  Fig.  2).  r^Cx)  is 
normally  distributed  from  our  second  hypothesis  and,  since 

& is 


dx 


linearly  related  to  rx(x) , it  also  holds  a Gaussian  distri- 


" 'tl 
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button.  Hence,  the  vector 


r (x) 
x 


dr 

dx 


-(x) 


is  normally  distributed. 

Therefore,  if  p(a,0,x)  represents  the  probability  that: 


a < rx(x)  ^ a + da 


dr 

8 « < 8 + d9 


(4) 


we  have: 


p(a,0,x)  - ^ |m|  1 exp  (-  j a 0 H-1 


a 

0 


(5) 


where: 


KO.O) 


M - 


ft 

H2 


(0,0) 


(6) 


see  Ref.  1. 
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Relation  (5)  can  easily  be  reduced  to: 


p(a,0;x) 


-1/2 


exp 


(7) 


In  which  we  have  set  for  convenience: 


mil 

Ko 


4 <«.«> 

H 


v - Y(0,0) 

O 


(8) 


Now,  let  us  find  the  probability  that  rx(x)  has  a zero  between  x 
and  x + dx;  in  the  interval,  all  the  cross-sections  of  the  possible 
realizations  of  the  ice  surface  are  regarded  as  straight  lines. 


At  first,  we  will  derive  the  probability  for  rx(x)  t0  have  a 
zero  with  positive  slope.  For  this  to  happen,  a and  0 have  to  be 


x + dx 


related  by  the  inequality 
(clearly  shown  on  the  sketch 
beside) : 

-0dx  < a < 0 
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Then , the  probability  for  rx<x)  to  have  a zero  between  x and 
x + dx  with  a positive  0 is: 


d6  I p(a, 0,x)da 


And,  the  probability  density  for  this  happening  is: 


d(Ax) 


« 0 
/ d0  / p(a,0, 

0 -0Ax 


x)da 


Ax  - 0 


which  reduces  to: 


(-0)  p(-0Ax,0,x)  d0 


Ax  * 0 


or: 


e p(O,0,x)  d0 


(10) 
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By  Che  same  token,  Che  probability  density 
zero  between  x and  x + dx  with  a negative  slope 


for  r (x)  to  have  a 
x 

is: 


-I 


6 p(O,0,x)  d0 


(11) 


Adding 
density  for 


relations  (10)  and  (11),  we  obtain  ff(x),  probability 

r (x)  to  have  a zero  between  x and  x + dx: 
x 


n(x) 


^ |0|  p(O,0,x)  d0 
•00 


Using  relation  (5)  and  Integrating,  it  comes  out: 


ir(x) 


(12) 


(13) 


The  Integral  of  ir(x)  over  a unit  length  gives  the  number  of  expected 
zeros  per  unit  of  a x-direction  since  the  integral  from  -®  to 
obviously  gives  the  total  number  of  zeros. 

Therefore,  Xx,  the  average  distance  between  two  zeros  of  rx(x) , 
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SUB -APPENDIX  C-2 


The  goal  of  this  appendix  is  to  calculate  Cy  defined  as: 


<\.  " < u<<x  <*«.*>>  - <u(<xm«v>»  <ymav>)>  > (y 

V)  qiax  loflx  max  max 


where : 


u(x.y) 


x/2  y/2 

/ /h"("’ 


U)  du  dl) 


-x/2  - y/2 


Expanding  the  square  in  Eq.  (1),  we  obtain: 


% - <0  <<x  >,  <y  >)>  - <2u(<x  >,  <y  >)> 

Vi  max  max  max  max 


<U(<x  >,  <y  >)>  + <u(<x  >,  <y  >)> 

max  'max  max  * max 


2 2 
“ <v^  (<x  >»  <y  >)>  - <u(<x  >,  <y  >)> 

U max  max  max  max 
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v (5  n)  • <h  (x  + £,  y + n)h”U,y)> 
i * 


i. 


■I 
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From  relation  (2),  It  Is  clear  that: 


<0(<x  >,  <y  >'>1 2  - <S  >2  <h“(C,n)>2  - V2  (9) 

mav  f ^mflv  mav  ' ' ' * 


max  'max 


max 


We  also  notice  that: 


<x  >//2  <y  >1/2 
max  ’ max 


<S  >2  <h"(f;,n)>2  - 16 

max 


/ 


<x  > 
max 


/2 


- - e 


<y  > 
max 


^2 


- n 


d£  dn 


Therefore,  from  Eqs.  (4),  (8)  and  (9),  we  end  up  with: 


, 


<x  >1/2  <y  >//2 
max  'max  _ 


max  'max 

— _ 

r c 

<x  > 

Cy  > 

/ / 

max  P 

max 

/ 

- £ 

n 

J J 

2 

2 

0 0 

L»  — 

L-  -J 

t- 


v (S.n> 


(10) 


[ 

i 


1 
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SUB-APPENDIX  C-3 


In  this  appendix,  ip  (t)  in  terms  of  tp( t)  will  be  derived. 


i|>  (t)  ■ <h  (x)h  (x  + r) > 

<Kt)  ■ <h(x)h(x  + t)>  (1) 


where  h (x)  is  defined  in  relation  (5) 

At  first , we  introduce  the  following  notations: 


h(x)  - a h (x)  ■ a 

h(x  + t)  ■ b h (x  + t)  ■ b (2) 

Let  D(a,b,T)be  the  first  joint  probability  density  of  the  couple 
(a,b).  D(a,b,l)  dadb  represents  here  the  percentage  of  times  that 
h(x)  will  lie  between  a and  a + da  while  h(x  + t)  lies  between  b 
and  b + db. 

Then,  from  our  second  hypothesis,  D(a,b,x)  is  Gaussian  with 
mean  value  0,  and  can  thus  be  written: 

D(a,b,i)  - ^ |m|_1  exp  (-  ^ t ||r||  M-1  ||r||  ) (3) 
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D(a,b,x) . 
It  Is: 


<Q> 


Q(a,b)  D(a,b,T)  dadb 


(8) 


Or,  when  identifying  Q with  the  product  a b 


Q 


a b 


when 

when 


a and  b < 0 
a or  b ^ 0 


(9) 


and,  substituting  for  D(a,b,x): 


a b 


*"(  t) 


2’<*02 


(a2  + b2)  ^ - 2ab\J> 

o 

2(\j)  2 - ty2) 

o 


Jdadb 


(10) 
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in  which  we  have  set  for  conveniences 


4,(0)  - * 


(ID 


Introducing  the  new  variables  x and  y such  that: 


a - x - y 

b - x + y 


(12) 


I 

1 

* 

I 


Equation  (10)  becomes,  noticing  that  the  Jacobian  of  the  trans 
formation  is  2: 


(13) 
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And,  making  use  of  the  error  function  defined  as: 


& 


erf  (x) 


x 

-[■-' 

/n  J 


dt 


(1A) 


we  end  up  with: 


*“(t>  - 


— 1 “ 

1/2 

r °° 
/ 

0J>o  - to 

4 

ir(i|<  + <|>) 
o 

J 

0 

w 

2 1 


(u  - 


exp 


( 2 ♦.  - *\ 
\ “ K + y 


y)  erf  (u) 


du  + 


*o  + » 

4/if  4> 


(15) 


We  notice  that  for  the  two-dimensional  auto-correlation  function 
vp“(^>n) , the  result  is  exactly  similar  to  Eq.  (15). 


’(5,n)  ■ <h"(x,y)h"(x  + C»  y + n)> 


KS,n)  ■ <h(x,y)h(x  + 5,  y + T<)> 


(16) 
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With  Che  notation: 


¥q  - ¥(0,0) 
¥ - ¥<5.n) 


we  have: 


¥"(5,n) 


(V  - f) 
o 

ir(Y  + 40 
o 


n 1/2 


- |)  erf  (u) 


f , * - *\ 


¥ + 4*1 


**«*  r 


r I«9l«.» 
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SUB-APPENDIX  C-4 


The  purpose  of  this  appendix  is  to  examine  the  behavior  of 

<|T(x)- 

From  relation  (15)  of  Sub-Appendix  C-3,  we  have: 


We  can  split  the  integral  into  two  parts: 


I. 


<***» 


11,1 7,  "tvyv'w  •rf'W‘*  ■•*'*?*.  ■ vr^jm 


After  some  algebra,  It  follows  that: 


- k 


n 


3 V2 
2 ‘ *#  - * 


exp  I - n 


♦o  - ♦ 
*0  + * 


2 ip 


Vi 


0Pq  + ip) 


16(i^0  -\p  ) 


erfcln 


»o  " ♦ 
*0  + ♦ 


M5) 


here,  erfc(z)  is  the  complementary  function  of  erf(z) 

1.  T - 0 

When  x ■ 0,  \p  - \p  and  therefore  relation  (5)  of  Sub- 
o 

Appendix  C~1  becomes: 
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having  n increase  to  infinity,  it  follows  from  (1),  (2),  (3)  and 
(6)  of  Sub-Appendix  C-l  that  the  only  remaining  non  zero  term  is: 
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2. 


f -►  00 

We  have  to  go  back  here  to  relation  (10)  of  Sub-Appendix 
C-3,  namely: 
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T -*■  °°  is  equivalent  to  have  \p  ■*  0 since,  at  infinity,  iJ>(t) 
is  zero. 

Hence,  integrating  the  above  relation  for  ip  ■ 0,  one  easily 
obtains: 
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which  is  not  surprising  since,  for  an  infinite  t,  h (x)  and 


h”(x  + x)  are  no  more  correlated  and  therefore: 


<h"(x)h“(x  + T)>  - <h“(x)>  <h"(x  + t)> 

■ <h'(x)>2  - 5} 

The  aspect  of  a typical  iT(x)  is  sketched  on  Fig.  8. 

3.  application 

For  numerical  computation  of  ij>  (x)  from  <Kx) , it  is  interest- 
ing to  notice  that: 

n - 3 in  relation  (3)  of  Sub-Appendix  C-l  is  highly  sufficient 
since: 
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implying  k3  * 1. 
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SUB-APPENDJX  C-5 


The  purpose  of  this  appendix  Is  to  give  some  details  on  the 
derivations  relative  to  the  mathematical  analysis  of  the  one-dimen- 
sional case.  Thus,  the  first  statistical  moment  <x  > of  the  pool 

max 

extent  and  its  variance  a will  be  determined  in  terms  of  a,  , the 

x h 

variance  of  the  roughness  height  and  i^(t)  the  auto-correlation  of 
the  latter. 

We  Introduce  the  notations  h+(x)  and  h (x)  such  that 
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and 
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h > 0 
h 4 0 
h < 0 

h * 0 (2) 


We  define  u(x) , amount  of  oil  per  unit  width  that  would  fill  the 
pockets  between  abscissae  0 and  x: 


o(x) 


MO  d£ 


(3) 


Preceding  page  blank 


209 


The  extent  of  the  pool,  x , for  a given  realization  of  the 

max 

ice  profile  is  such  that: 
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If  <x  > denotes  the  ensemble  average  of  x , we  can  always 

max  max 

write  Eq.  (A)  as: 
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since  h is  a continuous  function  of  x;  € corresponds  here  to  a 
point  of  the  profile  where: 
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As  was  done  for  the  two-dimensional  case  and  with  the  same  limitations, 
we  assume  that: 


h (O  - <h  (0> 


(7) 
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Equation  (5)  Is  then  written: 


h"«)  - <x  >)  <h"(0>  (8) 

max  max 


average  of  Eq.  (8)  and  noticing  that  this 
with  linear  operations: 
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Taking  the  ensemble 
kind  of  average  commutes 
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Then: 
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To  derive  the  variance  o of  x . defined  as: 

x max 


a - <(x  - <x  >)2> 

x max  max 
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we  come  back  to  Eq.  (8)  and  write  it  as: 
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Or: 
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where  a is  the  variance  of  the  volume  U (<x  >) : 

0 max 
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Also,  expanding  Eq.  (14): 
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From  Eq.  (3),  it  is  clear  that 
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where  ^ (t)  Is  the  auto-correlation  for  h . 
Making  the  following  change  of  variables: 


T - U - n| 

e - u + nl  (17) 


We  end  up  with: 


<V^(<x  >)> 

max 


(18) 


Noticing  that: 


(19) 


and  using  Eqs.  (13)  and  (15),  we  obtain: 


(20) 


if(T)  in  terms  of  ^(x)  is  determined  in  Sub-Appendix  C-3 
(Eq.  (15))  while  the  value  of  h~  has  already  been  determined  in 
terms  of  ah  in  Eq.  (22)  of  Appendix  C. 
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